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ABSTRACT 


A  method  haa  been  developed  for  computation  of  radiant  inter¬ 
change  within  an  electric  arc  where  self  absorption  is  important.  Moat 
or  all  radiation  in  the  vacuum  u.  v*  range  is  reabsorbed  before  it  gets  out 
of  the  arc.  Present  values  of  absorption  coefficients  are  within  an  order 
of  magnitude  of  predicting  radiation  losses  as  determined  by  direct 
experimental  measurement. 

An  arc  has  been  operated  in  air  at  100  atmospheres  sufficiently 
quiescent  to  measure  a  temperature  profile  at  100  amperes  but  great 
care  must  be  taken  to  start  symmetrically  and  avoid  disturbances.  At 
400  amperes  even  extreme  precautions  are  not  sufficient  to  forestall 
instabilities.  The  air  arc  constrictor  is  uncooled  and  is  designed  for 
quasi-steady  operation  for  1  to  5  milliseconds  depending  on  current  level. 


iii 


A  E  DC -TR -66-258 


NOMENCLATURE 


dA 

b 

C 


c 

c 


c 

p 

D 


d  , 

o 


dr  dz’  d3 


E 

f(x) 


g 

g 

H 


h 


h 

h 

I 

dthin 

J  ,  .  , 
thick 

k 


k 


(k  R) 


L 


Area  element 

dimension,  Figure  44 

constant,  Equation  2.  39 

radiation  constant  h/k  =  1.438  cm-°K 

B 

specific  heat  of  calorimeter  ring 
specific  heat  at  constant  pressure 

Relative  Cumulative  Spectral  Radiance  Equation  2,4 
Interpolation  Differences  Equations  A„  3-A.  6 
Voltage  gradient 

Transmission  integral  Equation  2*  16 
Calorimeter  Integral  Equation  Dt  2 
Transmission  function  Equation  2,  17 
dimension,  Figure  44 
dimensionless  enthalpy 
enthalpy 

calorimeter  half  width,  Figure  44 

t 

Plancks  constant 
Electric  Current 
Black  Body  Intensity 

Radiance,  Optically  thin,  Equation  3.  1 

Radiance,  Optically  thick,  Equation  3,  2 

number  of  data  points,  Equations  2,  36  and  2,  37  only 

linear  absorption  coefficient 

Thermal  conductivity  of  Air 

Boltzmann  Constant 

Thermal  Conductivity  of  Copper 

optical  depth 

length  of  transmission  path,  Figure  1 


v 


A  E  DC -TR -6  6*25  8 


Lo 

m 

m 

N 

N 

n 

P 

P 

r 

P 

(q/A) 

(q/z) 

R 

g 

R. 

jqm 

r 

s 

s 


s 

T 

t 

5 

V 


V 

w 


B,  n 


X. 

i 


x 

X 


y 

z 


Loschmidt  Number 

mass  of  calorimeter  ring 

parameter  in  Equation  2.  6 

Average  particle  density,-  Equation  C.  6 

number  of  Annuli  for  radiant  interchange  model 

wavenumber  (reciprocal  wavelength) 

2 

net  power  into  unit  volume  (a  E  -P  ) 

-  r 

net  radiated  power  per  unit  volume 
pressure 

heat  flux  per  unit  area 

heat  flux  per  unit  length  of  arc 

gas  constant 

sector  radii  Equation  2.  24 
radius  of  arc  or  arc  annulus 
dimension,  Figure  44 

number  of  sectors  per  quadrant,  radiant 
interchange  model,  Equation  2.  19 

parameter  in  Interpolation,  Equation  A.  3 

temperature 

time 

Time  of  arc  extinction 

volume 

hWk  T 

D 

Monochromatic  emissive  power  of  a  black  body 
Mol  fraction 

dummy  variable  in  transmission  integrals 

Temperature  or  enthalpy  profile  exponent, 
Equations  2.  35  and  5.  1 

Temperature  exponent,  Equation  C.  1 

axial  distance 


vi 


AEDC-TR -66-258 


a 

C 

a 

a 

a 

0 

y 

A 

€ 

c 


e 

x 

v 

i 

TT 

P 

a 


u> 


Greek  Letters 

Thermal  diffusivity  of  calorimeter 
Angle  defined  in  Equation  2.  22 
absorptivity 

Integration  limits,  calorimeter  integral, 

Figure  44 

Angle  from  normal  to  radiating  surface,  Figure  1 

increment,  not  necessarily  uniform 

emissivity 

ratio  of  radiant  flux  per  unit  length  reaching  calori¬ 
meter  to  total  radiant  flux  at  arc  boundary 

Angle  defined  on  Figure  1 

Wavelength 

F  requency 

Function  defined  by  Equation  2,  40 

3.  14. . .  . 

density 

Stefan  Boltzmann  radiation  constant 

Electrical  conductivity 

Standard  deviation 

transmissivity 

Angle  defined  on  Figure  2 

Solid  angle 


vii 


A  E  DC -T  R -6 6-25 8 


Subscripts  and  Superscripts 


a 

Zone  1  arc  boundary- 

C 

calorimeter 

CL 

centerline  value 

cond 

conduction 

e 

emitted 

i 

inner  boundary  of  an  annulus 

L 

Inner  boundary  of  outermost  annulus 

N 

Arc  core 

o 

reference  state  (1  atmosphere,  273°K) 

o 

outer  boundary  of  an  annulus 

P 

Plexiglas 

rad 

radiation 

t 

transmitted 

W 

wall 

Indices 

i 

general  summation  index 

i 

annulus  (counting  inward  from  the  outermost) 

j 

boundary  of  emitting  surface 

.  / 
j 

complement  of  j,  see  section  Zm  3 

m 

angle  increment 

n 

frequency 

q 

boundary  of  absorbing  surface 

vm 


A  E  DC  -TR-66-25B 


TABLE  OF  CONTENTS 

Page 

1.  INTRODUCTION . 1 

2.  ANALYSIS . 3 

2.  1  Model  and  Basic  Assumptions  .  •  .  . . . .  3 

2.2  Radiant  Interchange  Integral.  . .  4 

2.  3  Sector  Radii . . . • . . .  8 

2.4  Radiant  Interchange  Summation  Equations  ....... .  9 

2.  5  Two  Zone  Arc  Model  . . .  .  .  12 

2.  5.  1  Temperature  Profile . .  . .  12 

2.5.2  Thermal  Conduction . 13 

2.5.3  Electrical  Properties  . . *  15 

3.  AIR  ABSORPTION  COEFFICIENTS  AND  RADIANCE  .  17 

3.  1  Data  Sources  and  Ranges . .  »  * .  17 

3.2  Integral  Results  . . .  •  18 

4.  EXPERIMENTAL  APPARATUS .  22 

4.  1  Arc  Constrictor  Configurations.  . .  22 

4.  1.  1  Design  Principles . . . .  22 

4.1.2  The  Heat  Sink  Type  Cascade  Constrictor  •  •  .  •  . .  2  3 

4.  1.  3  The  Ablation  Type  Constrictor  . . .  24 

4.  2  High  Pressure  Test  Apparatus  •  •  •  .  • . .  . .  25 

4.  3  Power  Supply . 25 


IX 


AEDC-TR-66-258 


Page 

4.4  Instrumentation .  26 

4.4.  I  Electrical . • .  26 

4.4.2  High  Speed  Photography .  26 

4.4.3  Time  Resolved  Spectrograph  . . .  27 

4.4.4  Total  Radiation  Calorimeter .  29 

5.  METHODS  OF  MEASURING  TEMPERATURE  . 33 

5.  1  Continuum  Radiation  Intensity*  . .  33 

5.2  Arc  Temperature  From  Voltage  Gradient .  34 

6.  ARC  STABILITY .  35 

7.  RESULTS  AND  CONCLUSIONS .  37 

7.  1  Electrical  Gradients  and  Preliminary  Arc 

Temperature  Determination .  37 

7.2  Temperature  Measurements  From  Continuum  Intensities  ...  37 

7.2.  1  The  115  Ampere  Arc .  37 

7.2.2  The  400  Ampere  Arc . 38 

7.  3  Application  of  ARCRAD  III  to  Measured  Profile  .  38 

7.  4  Total  Radiant  Heat  Flux  Measurements .  39 

7.5  Conclusions .  41 

BIBLIOGRAPHY .  43 

APPENDICES 

A.  TABULATION  AND  INTERPOLATION  OF  TRANSMISSION 

INTEGRAL . 48 

B.  PROGRAM  ARCRAD  III .  54 

B.  1  Purpose . 54 


x 


AEOC-TR-66-258 


Page 


B.  2  Glossary  of  Terms .  54 

B,  3  Input .  64 

B, 4  FORTRAN  Listing,  ARCRAD  III  and  Subroutines .  67 

C.  ANNULUS  BOUNDARIES .  73 

C.  1  Weighted  Average  Temperature . 73 

C.2  Outermost  Annulus . .  •  •  .  74 

D.  TOTAL  RADIATION  CALORIMETER  MODIFICATION  TO 

RADIANT  INTERCHANGE  EQUATIONS .  76 

FIGURES .  79 


xi 


AEDC-TR-66-259 


1.  INTRODUCTION 

Planned  development  of  very  high  pressure  electric  arc  heaters  for  re¬ 
entry  simulation  facilities  has  resulted  in  a  need  for  reliable  calculation  and 
measurement  of  the  radiant  heat  transfer  from  arcs  operating  at  pressures 
of  the  order  of  100  atmos  ph  eres,  At  this  pressure  level  radiation  is  a 
critical  loss  mechanism  from  the  standpoints  both  of  arc  heater  efficiency 
and  container  wall  cooling  load.  Depending  on  frequency,  the  plasma  ranges 
from  mildly  self-absorbed  to  nearly  black  body,  but  Over  most  of  the  spectrum 
self-absorption  is  important  and  the  upper  bounds  on  radiant  intensity  fur¬ 
nished  by  the  optically  thin  and  black  body  approximations  are  completely 
inadequate.  Experimental  data  on  the  properties  of  plasmas  in  the  pressure 
range  around  100  ATM  is  scarce,  though  some  information  is  available 
from  the  recent  paper  of  Bohn,  Beth  and  Nedder  (Bo,  1 )  *\  the  well  known 
papers  of  Anderson  { A,  1)  and  Peters  (Pe.  1 )  and  a  few  more  (e.  g.  ,  Refs, 

Og,  1 ,  Mi.  1 ,  Bo.  1 ). 

The  objective  of  the  combined  analytical  and  experimental  program  de¬ 
scribed  in  this  report  was  the  determination  of  radiant  heat  flux  from  an 
axially  symmetric  air  arc  column  operating  at  high  pressure.  In  an  oper¬ 
ational  arc  heater  the  arc  column  is  likely  to  be  unconfined  or  only  partially 
confined  and  its  characteristics  are  likely  to  be  further  complicated  by  the 
presence  of  a  transverse  flow  and  perhaps  also  by  an  alternating  current 
power  supply,  but  restriction  to  a  steady,  axially  symmetric  arc  with  large 
length  to  diameter  ratio  and  no  flow  was  necessary  to  make  both  experiment 
and  analysis  manageable. 

An  analytical  model  for  calculation  of  radiant  flux  density  at  the  arc 
boundary  had  been  developed  under  a  previous  contract  and  this  served  as  the 
basis  for  a  more  advanced  model  with  which  local  net  radiant  interchange 
and  current  voltage  gradient  can  be  calculated.  The  model  considers  steady 
state  dependence  of  the  radiant  heat  flux  on  frequency,  temperature  and 
density,  including  the  effects  of  self-absorption,  by  breaking  up  the  arc  into 
a  series  of  concentric  constant  property  annuli  and  converting  volume  radiant 
emission  from  each  annulus  to  an  equivalent  emission  from  the  exterior  and 
interior  annulus  surfaces. 

The  experimental  arc  was  designed  to  provide  a  well  defined  arcing 
apparatus  suitable  for  detailed  studies  of  high  pressure  arc  plasmas,  based 

'^-References  are  listed  alphabetically  by  first  and  second  {occasionally  third 
or  fourth)  letter  of  the  senior  author1  s  last  name  and  a  number  to  distinguish 
more  than  one  publication  by  the  same  author. 
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in  part  on  previous  experience  with  Argon  filled  lamps  which  had  been 
operated  at  70  atmospheres  and  4000  amperes  for  one  millisecond*  It  was 
designed  for  short  duration  (~  700-800  ^s)  so  that  no  cooling  is  needed  de¬ 
spite  the  very  large  heat  flux,  thus  keeping  hardware  design  simple  and 
flexible.  Arc  power  supply  is  a  capacitor  bank  supplying  a  constant  current 
pulse  for  the  requisite  duration  by  means  of  an  RL C  network  switched  with 
spark  gaps* 

Analysis  and  experiment  represented  a  coordinated  effort  to  measure 
and  calculate  radiant  heat  flux  from  the  arc*  Radiant  heat  flux  both  leaving 
the  arc  and  reaching  a  calorimeter  were  computed  and  compared  to  calori¬ 
meter  measurements*  The  computations  were  based  on  a  directly  measured 
temperature  profile  and  on  analytic  profiles  constructed  with  some  guidance 
from  observations  of  electrical  properties* 

Though  computation  £.nd  experiment  both  depart  from  the  reality  of  an 
operational  arc  heater  both  have  given  useful  insights  into  problems  and 
possibilities  in  the  high  pressure  regime*  Keynote  of  the  experimental 
effort  has  been  the  struggle  for  operation  sufficiently  stable  to  obtain  useful 
spectrographic  data*  From  the  analysis  has  come  a  quantitative  indication 
of  the  remarkable  distinction  between  "low  frequency"  and  "high  frequency" 
energy  in  the  radiation  process. 
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2.  ANALYSIS 


Even  with  the  restrictions  of  steady  state,  axial  symmetry  and  large 
length  to  diameter  ratio,  radiant  emission  from  a  high  pressure  electric  arc 
is  an  extremely  complex  set  of  phenomena,  These  phenomena  are  character¬ 
ized  by  optical  depth,  the  dimensionless  product  of  an  absorption  coefficient 
(dependent  on  frequency ,  temperature  and  density)  with  an  orientation  dependent 
characteristic  length* 

The  objective  of  the  analytical  work  has  been  to  set  up  an  analysis  of  the 
radiant  interchange  within  the  arc  which: 

1)  exploits  the  data  storage  and  handling  capability  of 
a  large  digital  computer  to  take  all  these  variables 
into  account. 

2)  has  a  degree  of  refinement  adjustable  to  suit  require¬ 
ments  and  is  sufficiently  economical  of  computer  time 
that  many  cases  can  be  run. 

At  least  two  different  approaches  may  be  taken  to  the  problem.  These 
are  the  astrophysical;  in  which  the  dimensionless  optical  depth  is  taken  as 

the  primary  variable  in  a  set  of  non-linear  integral  or  integro-differential 
equations;  (Ko.  1),  and  the  engineering;  in  which  the  emission  from  a 
volume  is  converted  to  equivalent  radiation  at  a  surface  (Ja,  1).  Most  of 
the  work  in  the  former  category  has  been  done  for  the  plane  parallel,  or 
spherical  case  although  there  has  also  been  some  work  recently  on  cylindrical 
geometry.  Kesten  (Ke.  1)  has  applied  the  method  to  arcs  but  its  extreme  com¬ 
plexity  has  limited  it  to  calculation  at  only  a  few  specific  wave  numbers 
(Ke,  1),  DeSoto  (  De.  1)  has  calculated  radiant  transfer  in  a  pipe  flow  for  a 
limited  spectral  band. 

2.1  MODEL  AND  BASIC  ASSUMPTIONS 

The  engineering  approach  has  been  taken  on  this  project  and  the 
computational  model  is  a  substantially  revised  and  improved  version  of  one 
previously  reported  (Ma,  1,  Ma,  2),  The  radiating  arc  column  is  assumed 
to  be  cylindrically  symmetric,  with  a  length  to  diameter  ratio  sufficiently 
large  ( 1  0  or  more)  to  be  ‘'infinite".  The  arc  cylinder  is  divided  into  annuli, 
each  at  a  uniform  temperature  temperature  and  density.  Then,  taking  into 
account  the  effect  of  self-absorption  within  an  annulus,  volumetric  emission 
from  each  annulus  is  converted  to  an  equivalent  radiant  flux  from  the  annulus 
boundary  surfaces.  The  assumption  of  local  thermodynamic  equilibrium 
(LTE)  permits  the  application  of  Kirchoff's  law  relating  absorption,  trans¬ 
mission  and  emission.  Except  for  the  core,  each  annulus  has  an  interior, 
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concave,  surface  as  well  as  an  exterior,  convex,  surface  from  which  radiant 
emission  must  be  taken  into  account.  Radiant  interchange  among  all  annuli 
is  computed  by  considering  emission,  absorption  and  transmission  along  a 
representative  array  of  paths  from  all  the  emitting  surfaces  to  the  outside 
of  the  column. 

Because  self  absorption  is  significant,  the  controlling  parameter  is 
optical  depth,  the  dimensionless  product  of  absorption  coefficient  and  a 
characteristic  length.  Two  optical  depths  are  distinguished  here, 
l)  absorption  lengths  along  the  transmission  path  and  2)  emission  lengths 
within  each  annulus.  Since  absorption  coefficient  varies  with  temperature 
and  density  (each  a  function  of  radius),  absorption  and  emission  lengths 
are  computed  as  integrals  or  sumrra  tions.  For  example,  a  typical  absorption 
length  (Fig.  1)  is  given  by 

(kR)  =  f  kdr  =  sin  0  f  kdl  . 

Jo  Jo  (2.1) 

The  basic  geometrical  element  of  volume  used  for  calculation  of  emission 
and  absorption  is  the  infinitely  long  half- cylinder,  with  a  typical  element  of 
emitting  area,  dA  located  on  its  axis,  Fig.  1.  In  previous  work  (Ma.  1)  an 
average  value  of  emissivity  (frequency  dependent)  was  calculated  for  each 
annulus,  but  now  the  calculations  of  angle  dependent  emission  are  on  the  same 
basis  as  transmission. 


Neither  emission  nor  transmission  from  any  surface  element  actually 
takes  place  through  a  half  cylinder,  but  the  true  volume  shape  is  effectively 
approximated  by  several  cylindrical  sectors  of  different  radii,  Fig.  2.  Two 
typical  paths  are  shown  in  Fig.  2f  with  emission  lengths  shown  as  solid 
lines  and  absorption  lengths  shown  as  dashed  lines.  Calculations  of  total 
emission  and  absorption  in  each  wave  number  interval  requires  integration 
over  2  77  steradians,  accomplished  by  integrating  separately  over  angles  $ 
and  0(Fig.  1).  The  sector  approximation  decouples  the  angular  dependence 
of  absorption  length  or  emission  length  on  <£  within  any  sector. 

2.2  RADIANT  INTERCHANGE  INTEGRAL 

Solution  of  the  problem  of  radiant  interchange  among  the  annuli  proceeded 
as  follows.  Consider  a  surface  element  dA  of  an  annulus,  radiating  with  an 
average  intensity  I  within  a  solid  angle  dCCover  a  wavenumber  interval 
An  in  a  direction  making  an  angle  y  with  the  normal  to  dA.  This  equivalent 
surface  emission  arises  from  the  volume  within  the  annulus  contained  within 
the  projection  of  du  back  through  dA.  Some  of  the  radiated  energy  will  be 
absorbed  by  other  annuli  along  the  path  and  some  will  escape  the  arc  column 
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entirely.  Then  the  radiant  heat  flux  emitted  at  one  annulus  surface  and 
reaching  some  other  surface  inside  or  outside  the  arc  is  given  by 


d  q  =  I  £  T  cos  y  d  to  dA 


(2.  2) 


where  1^  is  black  body  intensity,  e 
T  is  transmissivity  along  the  path. 


is  equivalent  surface  emissivity  and 


Black  body  radiant  intensity  within  a  particular  wave  number  interval 
can  be  written  as 


r  =2JL 

B  17 


[  D(n)  -  D  (n  +  An)  ] 


(2.  3) 


where 


D  = 


In  WB.  „  dn  £ 

c 


W„  dn 
o  B,  n 


W_  dn 
B,  n 


(2.4) 


and  D  is  conveniently  expressed  for  computer  calculation  as  {Pi.  1), 


D  s  1  - 


15 

4 

77 


2  4 

v_  ,  _v_  _  _v _ 

8  60  5040 


V6  \ 

+  272160/ 


v  <  1 


(2.  5) 


TT 


-  mv 

- — —  (  (  (  mv+3)  mv+6)  mv-6) 

m 


v  >  1 


(2.  6) 
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h  V 

v  k  T 

B1  (2.7) 

Referring  to  Fig.  1*  the  solid  angle  subtended  along  a  typical  path  is 


dUJ  =  (R/L)  d  0  d  =  sin  0  d  6  d  $ 

also 


(2.  8) 


cos  y  -  sin  0  cos  $  (  2.  9) 

The  radius  of  the  cylindrical  emitting  surface,  Fig.  2,  is  r.  so,  by 
symmetry  ^ 


dA  =  2  7J  r,  dz 
J 


(2.  10) 


Transmissivity  along  the  path  length  L  is 


T  =  e 
n 


(2.  11) 


where  subscript  n  denotes  wavenumber  dependence,  so,  using  Eq.  (2.  1) 


-  (kR)t  /  sin  0 
n 


(2.  12) 


The  subscript  t  on  (kR)  refers  to  the  transmission  path.  Frequency 
dependence  is  understood  but  not  explicitly  indicated. 

Emissivity  at  the  radiating  surface,  dA,  is  computed  in  an  analagous 
way  except  that  the  absorption  length  is  replaced  by  an  emission  length,  the 
dashed  and  solid  lines  respectively  in  Figure.  2. 


€  =  1 
n 


-(kR)  /  sin  0 

e  e 


( 2.  1  3) 
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Eqe  (2.  2)  can  now  be  rewritten  in  terms  of  radiant  flux  per  unit  arc  column 
length  as 


[D(n)  -  D(n  +  4n)j 

1-  e  '(T^>e  1  3in6]  [e 
sin  0  d0  cos  <£d  0^  ^2  U  r,  J 


/  sin  0 


(2.  14) 


Absorption  and  transmission  lengths  are  assumed  constant  within  each  sector, 
A&,  so  that  Eq„  (2.  14)  can  be  integrated  over  $  and  6  separately,  with  a 
factor  of  4  arising  because  the  integration  is,  by  symmetry,  over  1  ^ 4  the 
solid  angle* 

(q/z)ln,A*=  Zn  Ij0T*  [«»>  -  D(n  +  An>]  [sin(*  +  -^-)  -sln(*-^)] 

1,12  e-[(kR)t  +  (kB)J  /  sin  e  sin2  e  d  el 

°  J  (2.  15) 

The  symbol  (q/z)  refers  to  radiant  flux  per  unit  length  of  arc.  Superscripts 
t,  e,  a  will  specify  transmitted,  emitted  and  absorbed  respectively. 


The  remaining  integrals  are  functions  of  optical  depth,  not  integrable 
analytically  but  identical  in  form 


f  (x)  = 


-x/sin  0 
e 


sin  0  d  0 


(2.  1 6) 


A  numerical  integration  has  been  performed  and  a  table  of  differences  pre¬ 
pared  so  tnat  the  function  can  be  evaluated  simply  for  any  x,  see  Appendix  A. 


Defining  a  transmission  function  in  terms  of  Eq.  {2.  16) 
g  =  f  [  <k"R)t  ]  -  f  [(  kR)t  +  <k“R)e] 


(2.  17) 
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Equation  {2.  15)  becomes 


(q/z) 


t 


21I  r.a  T 
J 


D(n  +  A  n)J  £sin(<£  +^pj-sin^ 


(2.  18) 


2.3  SECTOR  RADII 

Approximation  of  absorbing  or  emitting  volumes  as  sectors  of  cylinders 
(Fig.  2)  requires  calculation  of  the  radial  distance  from  each  annulus  boundary 
surface  to  every  annulus  boundary  surface  between  it  and  the  outer  surface  of 
the  arc.  The  distances,  measured  at  angles  0<  <£><7T/2  from  the  normal 
to  the  surface,  are  referred  to  as  sector  radii.  Surfaces  are  numbered 
through  the  arc  from  1  to  2  N,  where  N  is  the  total  number  of  annuli  includ¬ 
ing  the  core.  Sector  radii  are  then  assigned  a  triple  subscript  denoting: 

1}  the  number,  j,  of  the  originating  surface;  2)  the  number,  q,  of  the 
terminating  surface,  3)  an  index  m  on  angle  .  Examples  are  shown  in 
Fig*  3  for  the  case  N  =  3. 

It  should  be  noted  that: 


1}  The  number  j  is  never  less  than  q,  and  when  q  -  j,  R  =  0. 


2}  The  annulus  radii  r.  and  i* _  are  equal  .  They  will  be 

l  2N+1-1  7 

referred  to  as  complements  and  denoted  by  primes. 

(e.  g.  if  N  =  3,  r  =  r  /  and  r  =  r  /). 


3)  Physical  properties  and  radiant  emission  of  an  annulus  are 
referred  to  by  the  subscript  of  the  equivalent  emitting  surface,  so 
that  (q/z).  and  ( q / z) ^  refer  to  the  same  annulus  and  are 

eventually  combined,  (e.g.  if  N  =  3.  T  =  T  .  and  k  =  k  ) 

2  4  n,  1  n,  5 

4)  With  s  the  number  of  sectors  into  which  the  quadrant  0  <  4><  tt/2 
is  divided,  and  therefore  also  the  number  of  angles  at  which  sector 
radii  are  computed, 


and,  with  1  <  m  <  s 


A  $  -  IT  /  2s 
4>  =  (m  -4“) 

m  X  2  / 


*  +  A$ 

m  — 

cos  4  d  $  =  sin  (m  A  4-  )  -  sin  [(m-1)  A^  ] 

I  A4> 

$  -  — — 

m  2 


(2.  19) 

( 2.  20) 


(2.  21) 


8 


AEDC-TR.66-258 


Trigonometry  for  typical  sector  radii  is  shown  in  Fig.  4  for  the  three 
possible  relationships  between  j,  q  and  N.  For  convenience  in  calculating 
the  sector  radii,  an  angle  a  measured  from  the  outward  normal  of  the 
boundary  surface  is  defined  such  that 


$  =  a 

$  =  rr  -  a 


a  <  77/2 
a  >  77  /2 


{2.  22} 


For  either  Fig.  4a  or  4b  the  sector  radius  is  related  to  the  radii  of  the 
originating  and  terminating  boundaries  by 


/r.  coset  +  +  /r,  sina\^  =  r  i 

\  jqm  )/  \  jqm  /  q 

The  sector  radius  is  then  given  by 


R.  -  - r .  cos  3,  + 
jqn'  J 


J~2  27  ~~1.  2~ 

*  r .  cos  C.  j-  ( r  -  r .  ) 

'  )  q  j 


q  <  N 


{2.  23) 


R . 
jqm 


=  -r  cos  a,  -  V  r . 
j  T  J 


2  2  2  2 

cos  a+(r  -r.)  q>  N 

q  J 


(2.  24} 


with  the  condition  that 


r  >  r.  sin  3 

q  “  J 


(2.  25) 


If  condition  (2.  25)  is  not  met,  the  sector  radius  does  not  exist  for  that 
combination  of  j,  q  and  m.  (e.g.  ,  Fig.  3,  R^  _  does  not  exist.  This  means 
that  the  radius  of  surface  q  is  too  small  for  the  path  originating  at  surface 
j  at  angle  to  intersect.  When  this  situation  occurs,  appropriate  deletions 

must  be  made  from  the  summation  equations  to  follow. 


2.4  RADIANT  INTERCHANGE  SUMMATION  EQUATIONS 

In  terms  of  sector  radii  and  the  appropriate  (frequency  dependent)  ab¬ 
sorption  coefficients,  the  transmission  functions,  Eq.  (2.  17),  become 


g. 

&jqm 


=  f  f(k  R).  1  -  f[(kR).  -  (k  R)., 

L  jqmJ  L  jq™  j  q™. 


(2.  26) 


where 
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—  y 

(k  R).  =r  4^  k.  (R..  -  R.  .  ,  ) 

jqm  i=q  1  31m  j,  1+ 1 ,  m 


{2.  27) 


Emission  lengths,  (k  R }j  /  are  computed  using  the  shortest  existing 

complementary  sector  radius,  either  R ,  or  R./  .  .  For  example, 

J  .  J  - 1 ,  m  J  jm 

referring  to  Fig.  3  emission  length  for  radiation  along  sector  radius  R^.^ 
is  computed  using  R  (j-5  ,  j  =2,  j  -1  =  1,  m  =  2). 

C*  1 

Energy  emitted  by  an  annulus  is  just  the  special  case  of  zero  trans¬ 
mission  length,  for  which  Eq.  (2.  26)  reduces  to 


g.  .  =  1  -  £  f  {k  R).,  1 

J  J  ni  L  ]  qmj 


(  2.  28) 


and,  within  a  particular  frequency  interval,  radiant  emission  per  unit  length 
from  each  surface,  j,  is 

(q/z)J  =  2  TT  a  Tj  [  D(n)  -  D(n  I'An)] 


0 

53  (  sin  (mA^  )  -  sin  [(m-1)  A  $  ]  }  [  g.  .  ] 

m  =  l  j  J  m 


(2.  29) 


Total  emission  from  any  annulus  is  found  by  adding  together  results 
for  the  j  and  2N-j  surfaces  (the  "exterior"  and  "interior"  surfaces  of 
the  annulus) 


<q/z)%  E  [(q/zf,.  +  (q/z>j_2N_i]  [o(n)  -  D(n  +  A  n)J 


l  <  i  <  N- 1 


{q/z)G  =  E(q/z)6  [  D  (n)  -  D(n  +  A  n)  ] 

n  n  j  =  N 


(2.  30a) 


(2.  30b) 


Frequency  summation  is  over  those  wavenumbers,  indexed  by  n,  which 
are  included  within  the  limits  .  2  <  h  y/k^T.  <15* 

These  expressions  for  radiant  emission  per  unit  length  include  self¬ 
absorption  within  the  emitting  annulus.  In  the  limiting  case  where  an 
individual  annulus  reaches  the  back  body  limit  within  a  particular  frequency 
interval,  the  function  g.  .  becomes  unity  in  that  interval  for  all  m  and 
EqP  (2.  29)  reduces  to  ^  the  back  body  emission  of  a  cylindrical  surface 

of  circumference  2  77  r.* 

1 
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Radiant  energy  emitted  at  one  annulus  and  absorbed  by  another  annulus 
is  computed  as  the  difference  between  energy  transmitted  to  and  through  the 
annulus.  Absorption  between  surface  q  and  surface  q  4-  1  of  emission  from 
all  surfaces  j  is  thus 


2N-1 

D(n)  -  D(n  +  An)].Z  .  T4  r. 

j=q+i  j  j 


S 

£ 

m 


-g. 

q+l,m  jqmy 


[  sin  m  -  sin  (m-  l)  A  $  ] 


(2.  31) 


Total  absorption,  by  any  annulus,  of  emission  from  all  other  annuli, 
analagous  to  Eqs,  (2.30),  is 


(q/z)a  =  E[(q/z)  .  +  (q/ z)*  -  .  3  [D(n)  -  D(n  +  An)] 

*  ”  q=i  q-  di\  -l 


l  n 


(2.  32a) 


1  <  i  <  N-  1 


(q/z)  *  =  £  (q/z)a  [  D  {n)  -  D  (n  +  A  n)  ] 
In  n  q=I\ 


(2,  3  2h) 


The  net  radiation  from  each  annulus  can  be  positive  or  negative  and  is 
given  by 


(q/z)radr  (q/z)e  -  <q/z)a 

ill  (2.  33) 

and  finally,  the  radiant  flux  per  unit  length  which  reaches  the  outside  surface 
of  the  arc  is  given  by 


(q/z) 


rad 


N 

=  £  (q/z) 

i=l 


net 


(2.  34) 


The  above  procedure  is  a  massive  sorting  and  summing  job  ideally 
suited  for  a  digital  computer.  The  fact  that  not  all  transmission  paths  inter¬ 
sect  all  annuli  results  in  the  absence  of  some  terms  in  the  above  equations. 
These  omissions  have  not  been  indicated  explicitly  but  are  taken  care  of  in 
the  scheme  for  indexing  the  summation  loops  of  the  computer  program  called 
ARCRAD  III.  Details  of  this  program,  including  a  Fortran  listing,  will  be 
found  in  Appendix  B, 
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2.5  TWO  ZONE  ARC  MODEL 

Application  of  the  foregoing  radiant  interchange  equations  to  a  specific 
arc  requires  additional  information  and  boundary  conditions*  If  the  temper¬ 
ature  profile  is  known  by  direct  measurement,  then  ARCRAD  III  can  be  used 
to  compute  radiant  interchange,  and  existing  information  on  transport 
properties  employed  to  calculate  voltage  gradients,  current  and  power 
dissipation  in  the  arc,  for  comparison  with  measurements  of  these  parameters 
If  the  temperature  profile  is  unknown,  then  one  might  hope  to  establish  it 
from  known  boundary  conditions  supplemented  by  measurements  of  electrical 
char  acte  r  is  tics . 

Both  radiant  intensity  and  electrical  conductivity  in  air  change  very 
rapidly  with  temperature,  so  that  the  lower  temperature,  outer  part  of,  the 
arc  contributes  relatively  little,  either  to  power  dissipation  or  to  radiant 
emission.  Also,  the  accurate  determination  of  temperature  profiles  in  the 
outer  region  is  very  difficult*  For  these  reasons  the  arc  is  divided  into  two 
zones  -  interior  {Zone  I)  and  exterior  (Zone  II). 

In  Zone  I,  once  the  temperature  profile  has  been  established,  as  dis¬ 
cussed  below,  the  annulus  boundaries  are  so  determined  that  the  weighted 
average  temperature  of  each  annulus  is  a  multiple  of  1000  K,  This  pro¬ 
cedure  results  in  a  representative  but  tractable  number  of  annuli  and 
eliminates  the  need  for  temperature  interpolation  of  absorption  coefficent 
data.  It  is  further  discussed  in  Appendix  C. 

In  Zone  II,  the  profile  is  assumed  linear  and  the  boundaries  of  annulus 
i  set  at  radii  corresponding  to  T.  -  500  K,  where  T,  is  a  multiple  of  1000°K 
greater  than  or  equal  to  3000  K.  XThe  outermost  annulus  includes  the 
temperature  range  from  2500  K  down  to  wall  temperature.  Molecular  ab¬ 
sorption  of  vacuum  uv  radiation  is  the  only  significant  radiation  effect  below 
2500  K  and  calculation  of  absorption  coefficients  for  this  annulus  is  a  special 
case,  also  treated  in  Appendix  C,  Inclusion  of  Zone  II  in  the  radiant  emission 
absorption  calculation  changes  the  Zone  I  results  by  less  than  1%.  However, 
the  vacuum  u.v.  absorption  in  Zone  II  is  strong  enough  in  most  cases  to 
prevent  any  of  this  radiation  from  reaching  the  wall  (see  Section  7.  3), 

2.5.1  Temperature  Profile 


An  analytic  expression  for  temperature  profile  in  which  the  fullr  :ss  of 
the  profile  is  governed  by  a  single  parameter  is 


T 


CL 


T 


a 


(2.  35) 
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For  any  x  greater  than  1.0  this  equation  satisfies  the  basic  condition 
of  zero  slope  at  the  centerline  and  generates  a  family  of  curves  ranging 
from  essentially  linear  variation  with  radius  (  x  =  1),  to  uniform  temperature 
{x  «  ). 


The  value  of  x  which  most  closely  maintains  voltage  gradient  constant 
in  the  radial  direction  was  the  criterion  which  we  used  to  compute  profile 
shape.  However,  b  oth  radiant  emission  and  thermal  conduction  change 
very  rapidly  with  temperature,  and  computed  voltage  gradient  is  very 
sensitive  to  these,  so  no  clear  cut  values  of  x  could  be  established  without 
some  guidance  from  experimental  measurements.  Experimental  temperature 
data  were  fitted  to  Eq,  {2.  35)  using  a  standard  least  squares  technique  (Sb.  1)* 
Application  of  the  least  squares  criterion  to  Eq.  (2.  35)  yields 

k 


2 

(r.  >  0)  (2.  36) 


SD 


k 


E 

i=l 


x  In 


r. 

i 

r 

a 


( 2.  37) 


where  k  is  the  number  of  data  points  (r.,  T.),  CT  is  the  standard 

1  l  oJJ 

deviation  and  T  was  chosen  to  minimize  (T-—  . 

a  SD 


2-5-2  Thermal  Conduction 


In  the  inner  zone  of  a  high  pressure  arc,  where  practically  all  the 
electrical  dissipation  takes  place,  energy  loss  is  predominantly  by  radiation, 
so  a  relatively  simple  approximation  was  used  to  estimate  thermal  conduction. 
In  the  outer  zone,  radial  conduction  and  absorption  of  vacuum  u,v,  radiation 
from  Zone  I  are  assumed  to  be  the  dominant  modes  of  energy  transfer. 


The  arc  is  assumed  to  be  steady  and  cylindrically  symmetric  with  no 
flow  so  that  the  Heller-Elenbaas  equation  (Sk.  1)  is  satisfied.  This  equation 
can  be  written  in  terms  of  enthalpy  gradient  as 


l 

r 


fl)  — 


(2.  38) 


where  P  is  net  power  into  a  unit  volume  due  to  the  combined  effects  of 
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electrical  dissipation  and  radiation.  Note  P  is  positive  when  radiant 
absorption  exceeds  emission  even  if  electrical  dissipation  is  negligible. 


Weber  (We.  1,  We.  2)  has  correlated  the  ratio  /  c  as  a  function 
of  enthalpy  (Fig,  5)  in  the  form  ” 


k  /  c  =  C  h 
A  p 


0.  5 


(2.  39) 

He  has  also  pointed  out  that  by  defining  a  function 
£  =  Ch1'5 

(2.40) 

one  can  simplify  the  first  term  of  Eq.  (2.  38).  Applying  this  transformation 


I  A 

r  dr 


+  1.  5  P  =  0 


(2.41) 


In  Zone  II  the  boundary  conditions  are 


r  =  r  h  2=  0 

w 


r  =  r  h  =  h 

a  a 

When  P  is  constant,  Eq,  (2.41)  yields  an  enthalpy  profile  for  Zone  II 


and  an  enthalpy  gradient  at  the  zone  boundary,  r  = 


(2.  43) 


A  positive  value  of  P  in  Zone  II  tends  to  reduce  the  enthalpy  gradient 
at  the  zone  boundary,  so,  as  a  first  approximation,  P  was  set  equal  to 
zero.  Conduction  heat  flux  at  r  =  is  then  given  by 
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(q/z) 


cond 


r  -r 


4  it  C  h 


1.  5 


3  In  (r  /  r  ) 
w  a 


(2«  44) 


Because  the  conduction  term  is  relatively  small  in  Zone  1,  and  by 
symmetry  must  vanish  at  r  =  0,  and  also  because  transport  properties  in 
the  high  temperature  zone  are  ill-defined,  heat  loss  from  the  arc  by  con¬ 
duction  in  Zone  I  was  assumed  to  be  proportional  to  (r/r  rather  than 
make  use  of  the  derivative  of  an  empirical  temperature  profile. For  any 
annulus,  i 


(q/z) 


cond 


(q/z) 


cond 
r  =  r 


&  ■ « 


(2.45) 


2.5*3  Electrical  Properties 

For  any  annulus,  i,  the  local  energy  balance  is 
rad  /  \  cond 


EZ  G.  = 


.  ■  (C  ■  w 


(2.46) 


where  (q/zK ra<^  is  given  by  Eq.  (2.  31)  and  (q/z)C°n<^  is  given  by  Eq,  (2.45). 
Conductance  G.  of  an  annulus  of  unit  length  is  defined  as 


G.  =  2  m  f  O'  rdr 
i  J  r.  e 

J  i-1 


(2,47) 


th 


where  r.  is  the  outside  radius  of  the  i  annulus  counting  inward  from  the 
arc  boundary.  Local  electrical  conductivity  cr  is  computed  from  Weber's 
correlation  (We,  1),  Fig.  6. 


CTj  [  mho/cm]  =  7.  35  x  1C  ^  H  ^ 

0,  22  0.  74 

ri 


0^  [  mho/ 

cm]  = 

.  338  P 

if  cr  <  a_ 

then  u 

=  a. 

1  2 

e 

1 

if 

then  a 

e 

(2.43) 
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8  2  2 

where  H  =  h/R  T  ,  R  T  =  7.874  x  10  cm  /sec  .  p  is  measured  in 

go  go 

atmospheres  and  local  values  of  enthalpy  are  obtained  from  the  temperature 
profile  (not  the  annulus  average  temperature  used  for  computing  radiant 
interchange). 

Equation  (2.46)  was  used  to  calculate  local  voltage  gradient.  Average 
gradient  for  Zone  I  is 


avg 


siH I  _  Nk* 

1  £  G  ’ 

i  i 


+  (q/z) 


cond 

r-r 


(2.49) 


where 


G  =  2 


4 


0*  r  d  r 
o  e 


and  the  weighted  rms  deviation  of  local  values  from  E  is 

avg 


SD,  E 


J?(e-E  f  G. 
_  J  1  1  avg  /  l 
®  r. 


Total  current  carried  by  the  arc  is  simply 


(2.50) 


(2.51) 


I  =  G  E 


avg 


(2.  52) 


16 


A  E  DC-TR-66-258 


3.  AIR  ABSORPTION  COEFFICIENTS  AND  RADIANCE 
3.1  DATA  SOURCES  AND  RANGES 

Equilibrium  radiant  emission  and  absorption  in  an  air  plasma  is  a 
complex  set  of  phenomena  including; 

1.  Free -free  continuum  radiation  (Br emsstrahlung) 

2,  Free-bound  deionization  continua 

3,  Molecular  bands 

4.  Molecular  absorption  of  vacuum  ultraviolet  radiation 

5f  Lines 

Items  lf  2P  and  3  were  included  in  the  calculations  of  Nardone,  Breene, 
Zeldin  and  Riethof  (Na  ,  1)  which  were  the  source  of  absorption  coefficients 
used  in  computations  with  ARCRAD  I  (Ma.  1). 

As  part  of  this  project  a  revised  tabulation  of  absorption  coefficients 
was  undertaken.  The  tabulation  is  a  composite  of  information  from  several 
sources  as  shown  on  Fig.  7.  Tjhe  Nardono,  et.  al.  results  were  retained 
at  wavenumbers  to  59,  000  cm  and  be^loy^  ^nce  they  are  most  complete; 
but  were  multiplied  by  the  factor  ( 1  —  e  B  to  correct  for  stimulated 

emission  (Sh.  l). 

From  59,  000  cm  ^  shown  in  Fig,  7,  the  recent  results  of  Sherman 
and  Kulander  (Sh.  1)  for  Free -Bound  radiation  were  used  together  with 
Kramer's  formula  (U.  1)  for  free -free  radiation. 

Sherman  and  Kulander  used  the  more  sophisticated  "quantum  defect" 
method  of  Burgess  and  Seaton  (Bg.  1)  rather  than  Breene's  hydrogenic  model 
to  calculate  radiation  due  to  radiative  recombination  of  singly  ionized  nitrogen 
and  oxygen  atoms  with  electrons,  and  these  are  the  dominant  contributors 
at  high  frequencies, 

A  quantum  mechanical  correction  to  Kramer’s  formula,  the  Gaunt 
factor,  was  assumed  to  be  unity  because:  l)  Nardone,  et.  al.  ,  used  unity 
in  their  calculation  of  Free-Free  radiation  at  low  frequency.  2)  The  number 
is  of  order  one  and  some  disagreement  exists  as  to  whether  it  is  greater  or 
less  than  one,  (No.  1)  (Ka.  1). 

Molecular  absorption  in  the  vacuum  ultraviolet  was  accounted  for  by 
averaging  the  experimental  data  for  molecular  absorption  collected  by 
Schultze,  et.  al.  (Sc.^l).  Molecular  absorption  takes  place  at  temperatures 
too  low  (  T  ^  10,  000  K)  for  equilibrium  emission  to  be  significant  at  high 
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frequencies,  but  the  absorption  process  is  an  important  one  in  that  vacuum 
ultraviolet  radiation  is  blocked  from  reaching  any  solid  boundary. 

The  effect  of  line  radiation  has  not  been  accounted  for  in  the  tabulation 
of  absorption  coefficients. 

Absorption  coefficients  on  the  inpi^t  data  tape  are  tabulated  at  75  waye- 
numbers,  equally  spaced  to  13,000  cm  and,  from  there  to  200,000  cm  , 
at  spacings  chosen  by  Sherman  and  Kulander  to  be  close  together  where 
absorption  edges  cause  sharp  changes  in  k  and  widely  separated  where 
changes  in  k  are  gradual.  Values  are  averages  over  the  interval  between 
the  tabulated  wavenumber  and  the  next  higher  tabulated  wavenumber. 

The  data  of  Nardone  were  available  at  decade  intervals  of  log  density 
ratio,  and  that  of  Sherman  and  Kulander  on  a  per  particle  basis.  Tabulation 
for  this  project  was  by  both  density  and  pressure  using  equilibrium^ thermo¬ 
dynamic  data  of  Browne  (B.  1).  Density  range  was  from  10  to  1  0  p/  p 

o 

and  pressures  tabulated  were  1,  5,  10,  30,  50,  100,  150,  200  and  300  at¬ 
mospheres.  Temperature  range  was  from  3000  K  to  20,  000°K  at  1000° 
intervals  with  an  additional  tabulation,  at  273  K,  of  molecular  absorption 
data  in  the  vacuum  ultraviolet. 

Fig.  8  is  a  sample  plot  of  k  for  several  temperatures  with  p  =  100 
atmospheres . 


3.2  INTEGRAL  RESULTS 


Condensation  of  the  enormous  amount  of  data  represented  by  frequency, 
temperature  and  pressure  or  density  dependent  absorption  coefficients  is 
highly  desirable.  This  can  be  done  by  summing  over  all  frequency  intervals, 
but  such  summations  must  be  treated  with  care  because  any  effort  to  include 
self-absorption  must  inevitably  be  associated  with  a  specific  geometry. 


Several  such  summations  have  been  performed  with  the  data  assembled 
for  this  project,  as  follows: 


1)  Radiance  from  an  infinitesimal  volume  element  (watts/cm  -ster); 


thin 


g  T 

71 


k  [  D(n)  -  D{n  +  An)  ] 
n  n 


(3.  1) 


This  form  is  independent  of  geometry  but  is  not  black  body  limited  and 
is  thus  applicable  only  when  the  gas  is  optically  thin  at  all  frequencies. 

Figure  9  is  a  plot  of  summation  Eq.  (3.1),  for  p  =  200  atm,  T  =  1  5,  000°K, 
and  gives  an  indication  of  the  relative  importance  of  any  wavenumber  interval. 
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TABLE  3.1 

COMPARISON  OF  RADIANCE  VALUES 


This  report,  composite 

Brsene  and  Nardone’’ 

(Na.  1} 

Yos  ^ 

(Yo.  1) 

Biberman*** 

(Bi.  1) 

J  ,  . 
thin 

3 

(kw/  cm  - s tor) 

J 

thick 

Lf  -  1  cm 

2 

( kw  /cm  -ster) 

J  ,  . 

thm 

(kw  /  c  m^-s  ter 

J 

thick 

L  =  1  cm 

2 

(kw /cm  -ster) 

\hin 

(kw/cm^-ster) 

^thick 

L  =  1  cm 

2 

(kw/cm  -ster) 

Pres  sure 

(atm) 

Temperature 

(°K) 

12000 

.  100 

.  080 

.  202 

.  140 

.  06 

.49 

1 

14000 

.  495 

.  379 

1.  26 

.  743 

.  35 

2.  22 

17000 

.  924 

*  858 

- 

- 

.  66 

4.  85 

12000 

1. 22 

.  625 

2.  47 

.  771 

.  86 

1.76  ! 

10 

14000 

8.  81 

3.  25 

21.  3 

4.  12 

7.  3 

6.  45 

17000 

35.  13 

13.  96 

- 

- 

35. 

28.  8 

12000 

1 

18.2 

6.  37 

27.  1 

6.  19 

- 

- 

100 

14000 

109.  2 

22.  9 

249. 

26.  4 

- 

45,  0 

17000 

656. 

78.  9 

- 

- 

- 

90.  7 

-(hiVkgT) 

*  not  corrected  for  stimulated  emission  factor  (1-e  ) 

**  from  curves  of  radiated  power 

from  tabulation  of  equivalent  surface  emissivities 
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Yos  (Yo,  1)  presents  the  results  of  his  calculations  as  total  radiated  power 
per  unit  volume,  i.  e.  ,  4  77  Table  3.  1. 

2 

2)  Radiance  from  a  slab  of  unit  thickness  {watts/cm  -ster,  L  =  1cm); 


thick 


4  /  -k  LA 

=  -?  T  —  2~>  (l.e  n  )[  D  (n)  -  D  (n  +  An)  ] 

77  n  \  / 


(3.2) 


Radiance  is  tied  to  the  doubly  infinite  slab  geometry  and  reduces  to  black 
body  intensity,  CT  T^/77,  in  the  limit  when  k L  >  >  1  at  all  frequencies.  As 
long  as  kL  >  >  1  the  magnitude  of  L  is  irrelevant  and  the  units  of  J^ick 
are  watts/cm^-ster).  When  kL  <  <  1  at  all  frequencies,  radiated  power 
is  proportional  to  slab  thickness,  L,  and  J 


thick 


is  identical  to 


J  -  in 

thin 


magnitude  and  units.  Nardone,  et.  al.  (Na,  1)  present  their  integrated 
results  as  radiance  from  a  slab  and  this  form  can  also  be  related  to  results 
of  Biberman,  et.  al.  (Bi.  1),  who  tabulate  equivalent  surface  emissivities 
for  a  slab,  Table  3.  1. 

3)  Radiant  Flux  density  at  the  surface  of  an  optically  thin,  uniform 
temperature,  infinite  cylinder  (watts /cm^); 


/  /  a  i  .  t  /  Volume  \  ^ 

(q/A)thin  ”  4  V  Jthin  (surface/  "  *  T a.  Jthin 


(3.  3) 


This  is  simply  the  application  of  J  ^  to  the  infinite  cylinder  geometry. 

4)  Radiant  flux  density  at  the  surface  of  an  optically  thick  uniform 
temperature  infinite  cylinder  (watts/cm  ); 


(q/A)thick  =  CT  T 


.  /  -1.9  k  r  \ 

na)[D 


].  9  k 

(n)  -D  (n  +  An)  ] 

(3.4) 

Here  the  mean  beam  length  concept  (Ec.  1)  has  been  applied  to  calculate  the 
equivalent  surface  emissivities  of  the  cylindrical  volume  and  thus  acccunt  for 
self-absorption.  Figure  10  shows  curves  of  (q/A)^.  ^  1 >  10,  100,  and  300 

atmospheres  for  ,  05  £  t  ^  1.0  (cm)  and  7000  ^  20,  000°  K  . 

In  Fig.  11  the  ratios  (q/A)^.^/  (<l/A)thick  and  (q/A)  /  O  T4  are 

plotted  as  a  function  of  pressure,  for  a  temperature -radius  combination 
typical  of  experimental  observations  at  100  atmospheres,  to  illustrate  the 
significance  of  self-absorption  in  arcs.  The  small  arc  radius  encountered 
at  100  atmospheres  tends  to  reduce  optical  depth.  Thus  self-absorption, 
while  significant  is  not  overwhelming,  as  shown  by  the  ratio  (q/A)  ,  .  / 


(q/A ) 


thick 


when  radiation  at  wavenumbers  up  to  59,  000  cm"1  is  considered 
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However,  when  vacuum  ultra-violet  radiation  is  included,  self  absorption 
is  very  important. 

Now  let  us  compare  (q/A)  to  the  black  body  limit.  As  pressure 

increases,  the  vacuum  u.v.  appears  less  important  because,  having  reached 
black  body  maximum  at  low  pressure,  its  relative  contribution  decreases  as 
low  frequency  radiation  becomes  stronger. 
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4,  EXPERIMENTAL  APPARATUS 


Our  experimental  efforts  to  operate  and  make  useful  measurements 
on  high  pressure  electric  arcs  operating  in  air  were  directed  toward 
achieving  steady  state  arcs  while  depending  on  the  transient  characteristics 
of  a  capacitor  power  supply  and  an  uncooled  constrictor*  This  permitted 
simplicity  and  flexibility  in  apparatus  design  but  required  fast  response 
instrumentation,  both  to  ascertain  if,  when  and  for  how  long  a  steady  arc 
was  established,  and  to  make  useful  measurements  within  the  steady  period* 
Figures  12  and  13  are  photographs  of  the  setup. 


Figure  14  is  a  schematic  which  shows  the  principles  of  operation. 
Time  resolving  element  for  the  temperature  measurements  is  a  rotating 
mirror  M  t  on  which  light  from  the  arc  is  focused  by  lens  L  ,  From  the 
mirror  the  light  falls  on  lens  and  is  then  focused  on  the  slit  of  a 
spectrograph.  Another  light  beam  from  the  arc  is  focused  by  lens  L  on 
slit  S  of  the  drum  camera,  which  provides  a  simultaneous  record  of  arc 
stability.  The  electrical  and  electronic  setup  for  the  generation  of  the 
high  current  pulse  and  for  the  proper  synchronization  of  the  measurements 
is  also  indicated  in  Fig.  14.  Details  of  all  the  experimental  equipment  and 
its  operation  are  discussed  below. 


4.1  ARC  CONSTRICTOR  CONFIGURATIONS 

4.1.1  Design  Principles 


Quantitative  information  as  to  plasma  properties  can  best  be  gained 
with  simple  arcing  devices,  specifically,  cylindrically  symmetric  ones 
without  flow.  A  particularly  successful  arc  constrictor,  in  which  cylindrical 
arcs  have  been  operated  for  long  periods  of  time,  is  the  so  called  cascade 
type  or  Maecker  arc.  (M.l,  M.  2).  This  constrictor  consists  essentially  of 
a  stack  of  mutually  insulated,  water-cooled,  copper  discs  with  a  central 
circular  hole  for  the  arc.  The  principal  limit  of  this  device  appears  to  be  ^ 
a  maximum  acceptable  steady  state  wall  heat  flux  of  the  order  of  1  0  KW/cm 
of  wall  surface,  (Co.  1,  Cu.  1,  Ga.  1).  Unfortunately  this  number  is  too 
small  to  operate  high  current  arcs  in  the  100  atm  range,  where  radiative 
losses  can  reach  values  which  are  10  or  even  100  times  as  high.  There 
was  hope,  however,  that  the  excellent  properties  of  the  Maecker  arc 
could  be  retained  by  operating  on  a  transient  rather  than  on  a  steady  state 
basis.  Then  the  copper  walls  could  be  used  as  a  heat  sink  and  the  arc 
switched  off  before  melting  occurred.  Figure  15  shows  the  result  of  a  rough 
calculation  of  the  maximum  permissible  operational  time  of  such  a  constrictor, 
using  conventional  heat  transfer  calculation  technique  (Ca.  1).  The  copper 
walls  can  be  expected  to  absorb  heat  flux  densities  of  100  KW/cm^  for  about 
1  millisecond  and  this  pair  of  parameters  was  attractive  enough  to  warrant  a 
further  study  of  such  a  device. 
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A  second,  quite  different,  solution  for  a  cylindrical  arc  constrictor 
which  can  be  operated  at  very  high  radial  heat  flux  densities,  is  the  ablation 
type  arc  (A,  1,  Bu.  1,  Fr.  1,  Og.  1).  This  is  a  constrictor  in  which  the  fluid 
or  solid  walls  evaporate  by  the  heat  from  the  arc.  If  proper  precautions 
are  taken,  the  arc  then  burns  in  the  ensuing  radially  symmetric  flow  field. 

From  a  theoretical  point  of  view  this  device  is  not  as  simple  as  the  cascade 
type  constrictor  when  the  latter  is  operated  with  no  flow.  However,  the  flow 
field  of  an  ablation  type  arc  is  well  defined  and  it  appears  that  a  number  of 
arc  properties,  especially  radiation,  can  be  studied  conveniently.  Moreover, 
a  recent  theoretical  treatment  of  this  arc  type  opens  the  road  further  for 
experimental  work  (Sm,  1).  We  have  built  both  types  of  constrictors  for  the 
purpose  of  investigating  high  current  air  arcs  in  the  100  atm  range,  but 
most  of  our  work  to  date  has  been  done  on  the  cascade  arc. 

4.1,2  The  Heat  Sink  Type  Cascade  Constrictor 

To  avoid  shorting  of  the  arc  at  the  high  electrical  gradients  encountered 
(about  200  v/cm)  the  copper  disc  elements  of  the  cascade  constrictor  must 
be  thin.  The  discs  used  are  20  mil  (1/2  inm)  thick  solid  copper  with  a 
central  hole  2,  0,  2,  5  or  5.  0  mm  diameter.  They  are  each  separated  by  a 
20  mil  air  gap  and  held  by  two  thin  ribs  of  Plexiglas  to  which  they  are  glued, 
Figs,  1 6  and  17.  During  the  assembly  operation  the  discs  and  electrodes  are 
aligned  by  a  steel  rod  through  the  constrictor  holes.  The  length  of  the  cascade 
type  arc  constrictor  is  7  cm,  thus  £/d  is  about  28.  As  can  be  seen  from 
Fig.  17,  the  copper  constrictor  is  radially  open,  permitting  a  fast  initial  ex¬ 
pansion  of  the  arc  plasma  after  ignition.  This  initial  flow  phenomenon  was  the 
slowest  transient  process  in  the  arc  and  in  a  radially  closed  constrictor  lasted 
several  milliseconds,  which  was  much  too  long  for  our  purpose. 

Elimination  of  water  cooling  has  several  important  advantages.  It  con¬ 
stitutes  a  remarkable  saving  in  design  and  operation,  permitting  economical 
experimentation  with  a  large  number  of  modifications  to  constrictor  diameter 
and  thickness  and  to  spacing  of  the  discs.  Also,  frequent  mishaps  such  as 
melting  of  the  copper  discs  by  faulty  arc  operation  can  be  borne  philosophically. 
Another  advantage  of  the  omission  of  water  cooling  is  the  convenient  access¬ 
ibility  of  the  arc  for  optical  observations.  This  helped  immeasurably  when 
arc  stability  was  being  studied  with  high  speed  photographic  techniques. 

Arc  electrodes  are  made  of  graphite.  They  are  held  centered  by  a 
thicker  copper  washer  at  a  distance  of  60  mils  from  the  first  copper  disc  of 
the  cascade.  For  additional  support,  and  to  aid  in  obtaining  a  good  electrical 
contact  between  the  graphite  electrodes  and  the  power  supply  leads,  the 
rather  light  constrictor  cage  was  mounted  in  a  strong  Plexiglas  holder  with 
two  copper  end  caps,  in  which  the  graphite  electrodes  were  held  smoothly  but 
tightly,  Figs.  16  and  17.  The  Plexiglas  holder  fits  into  the  pressure  vessel 
described  in  a  later  section. 
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To  improve  the  optical  quality  of  the  spectroscopic  measurements,  the 
two  central  elements  of  the  constrictor  were  modified  considerably  from 
those  described  above.  Fig,  17.  They  were  cut  to  a  rather  long,  nearly 
rectangular  shape  with  two  additional  tapers  designed  to  minimize  stray 
light  entering  the  optical  system  from  the  highly  reflecting  copper  walls. 

The  special  constrictor  elements  also  permitted  the  use  of  a  rather  large 
aperture  for  the  optics,  which  helped  in  obtaining  a  high  time  resolution 
with  the  photographic  spectroscopic  technique. 

Special  precautions  had  to  be  taken  with  the  copper  discs  to  prevent  arc 
contamination.  The  punchings  as  delivered  by  the  manufacturer  had  loose 
particles  and  small  greenish  crystals  on  the  surface  of  the  hole.  Also,  a 
burr  of  rather  loose  copper  material  was  apparently  unavoidable,  Figure  I  8. 
The  situation  grew  worse  when  the  die  used  for  punching  lost  it  sharpness. 
Loose  crystals  and  copper  particles  will  naturally  evaporate  much  more 
rapidly  when  exposed  to  the  arc  than  indicated  by  Figure  15.  Copper  vapor 
will  then  enter  the  arc  in  uncontrolled  amounts,  possibly  affecting  electrical 
conductivity  and  radiation,  and/or  giving  rise  to  local  unsymmetric  blowing 
which  may  cause  arc  instabilities.  To  avoid  these  effects,  all  copper  discs 
were  carefully  degreased  and  then  chemically  etched.  This  process  removed 
the  loose  particles  along  with  about  2  mil  of  the  copper  surface,  Fig.  18, 

Placing  of  the  fuse  wire  for  arc  ignition  in  the  copper  constrictor  was  an 
important  and  rather  time  consuming  operation.  To  minimize  contamination, 
tungsten  wires  as  thin  as  1  / 2  mil  (0.  00051’)  were  used.  For  maximum  arc 
stability  these  had  to  be  strung  straight,  and  accurately  centered,  in  the 
little  cascade.  This  was  achieved  by  replacing  the  two  graphite  electrodes 
having  a  central  hole  the  same  diameter  as  the  copper  discs,  which  were  used 
for  alignment;  by  two  others  with  a  tiny  central  hole.  Through  this  hole  the 
l/Z  mil  tungsten  wire  was  strung  and  held  taut  by  little  blobs  of  putty,  hung 
on  its  ends.  Then  the  wire  was  fixed  at  the  ends  of  the  electrodes  with  two 
more  little  blobs  of  putty. 

To  further  reduce  tungsten  contamination  of  the  arc,  1/4  mil  diameter 
wire  was  tried.  This,  however,  proved  to  be  such  a  tedious  task  that  it  was 
abandoned. 

A  good  indication  of  successful  (i.e.  reasonably  stable)  operation  of  the 
arc  was  the  state  of  the  copper  constrictor  after  arc  firing.  Symmetrically 
burning  arcs  did  not  melt  copper;  unsymmetric  arcs  burned  the  discs  badly, 
Fig.  19- 

4.1.3  The  Ablation  Type  Constrictor 


An  ablation  type  constrictor  was  also  built  for  some  exploratory  tests 
of  arc  stability.  To  check  on  stability,  complete  visibility  of  the  arc  is 
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advantageous  and  therefore  Plexiglas  was  used.  The  internal  diameter  of 
the  Plexiglas  tube  was  3  mm  and  its  length  3  cm*  In  the  form  shown  in 

Figs.  16  and  17,  without  the  use  of  special  windows  to  exclude  the  lens 
action  of  the  Plexiglas  tube,  the  assembly  of  the  constrictor  was  very  simple. 
Operation  was  also  simplified  by  the  radial  gas  flow  from  the  walls  which 
automatically  centered  and  stabilized  the  arc  so  that  location  of  the  ignition 
wire  was  not  critical. 


4.2  HIGH  PRESSURE  TEST  APPARATUS 

The  arc  constrictors  were  placed  in  a  strong  Plexiglas  holder,  and 
this  was  mounted  in  an  aluminum  pressure  vessel,  Figs.  20  and  21,  capable 
of  withstanding  at  least  5000  psi.  The  pressure  vessel  consists  essentially 
of  an  aluminum  olock,  into  which  several  openings  are  drilled.  First,  there 
is  a  large  central  hole  for  the  mounting  of  the  Plexiglas  holder  with  the  arc 
constrictor.  The  block  also  has  two  rectangular  holes  which  are  used  for 
arc  observation,  and  optical  measurements.  Finally,  there  are  two  more 
small  round  holes  for  auxiliary  purposes  such  as  locating  the  over-pressure 
relief  membrane.  Caps  for  all  holes  were  sealed  with  O-rings. 

Flow  of  electric  current  inside  the  pressure  chamber  is  schematically 
indicated  by  Fig.  22a.  The  vessel  i$  grounded  and  provides  a  return  path 
for  the  current.  This  helps  to  symmetrize  the  magnetic  field  lines  and  thus 
avoid  unsymmetric  magnetic  forces  on  the  arc.  To  further  reduce  stray 
light,  two  apertures  were  mounted  in  each  of  the  light  paths  before  the  large 
windows,  Fig.  22b. 

4.3  POWER  SUPPLY 

Power  is  supplied  to  the  arc  by  a  20  KV,  360  jJii  capacitor  bank.  A 
nearly  rectangular  current  pulse  is  produced  by  a  network  of  capacitors, 
inductors  and  resistors.  Essentially  three  parallel  circuits  are  used, 

Fig.  14;  a  weakly  damped  LC  circuit,  an  RC  circuit  and  a  strongly  damped 
LC  circuit  with  about  three  times  the  frequency  of  the  first  one.  The  cur¬ 
rents  of  the  individual  circuits  and  their  composite  are  shown  in  Fig.  2  3, 
which  also  illustrates  the  flatness  of  the  current  pulse.  Arc  resistance  has 
some  effect,  but  resulting  deformation  of  the  current  pulse  is  corrected  by 
small  adjustments  in  the  circuit  components. 

The  pulse  is  triggered  by  the  rotating  mirrpr  which  also  provides 
the  time  resolution  for  the  spectroscopic  measurements  so  a  rapid  and 
precisely  timed  start  of  the  current  must  be  ensured.  This  is  done  with  a 


25 


A  E  DC -T  R  -66-258 


tri-electrode  spark  gap,  G1  (Fig.  14)  fired  at  a  pre -determined  position 
of  the  rotating  mirror*  The  trigger  pulse  is  initiated  by  a  light  beam 
reflected  from  the  mirror  to  a  light  sensitive  transistor*  The  transistor 
in  turn  discharges  a  capacitor  through  the  primary  winding  of  a  pulse 
transformer,  delivering  a  steep  pulse  of  about  10  KV  peak  voltage. 

To  limit  the  dissipation  of  arc  energy  in  the  constrictor,  a  second 
spark  gap,G2  is  used  to  terminate  the  current  pulse*  This  gap  works  under 
rather  unusual  conditions  (Fig*  14)*  It  must  withstand  high  transient 
voltages  which  occur  both  at  the  time  of  ignition  of  gap  Gl,  and  during  the 
subsequent  explosion  of  the  tungsten  fuse  wire*  Later  the  gap  must  break 
down  with  only  a  limited  voltage  supplied  to  its  electrodes*  Therefore,  a 
special  trigger  gap  was  built.  Fig*  24,  with  a  separate  energy  supply  to  a 
fourth  electrode  which  generates  additional  plasma  for  the  breakdown  of 
the  gap*  The  design  of  the  4  electrode  gap  is  shown  schematically  in 
Figure  24*  Between  the  usual  insulated  central  trigger  electrode  E3,  and 
the  graphite  electrode  E2,  a  fourth  concentric  electrode  E4  is  mounted, 
slightly  recessed  into  the  Delrin  insulator*  The  four  electrode  gap  operates 
as  follows: 

When  the  initial  trigger  spark  flashes  from  E3  to  E2,  the  path 
between  electrode  E4  and  E2  also  breaks  down,  discharging  the  capacitor 
C*  This  discharge  creates  a  plasma  jet  by  evaporating  material  from  the 
narrow  slot  of  Delrin*  The  plasma  jet  is  directed  towards  El  and  assists 
greatly  in  breaking  down  the  path  between  El  and  E2,  even  when  only  a 
small  voltage  difference  exists  between  these  electrodes* 

4.4  INSTRUMENTATION 
4.4.1  Electrical 

The  arc  current  was  measured  with  a  coaxial  shunt  having  a  time 
constant  of  2*4  nanoseconds,  a  resistance  of  0.  01  ohm  and  an  1^  Rt 
value  of  500  joules*  Arc  voltages  were  measured  with  Tektronix  high 
voltage  probes  which  have  a  band  width  from  DC  to  10  M  Hz* 

All  electrical  measurements  were  monitored  with  Tektronix 
oscilloscopes  model  #551  or  #502A.  The  502A  instrument  was  used  for 
the  heat  flux  measurements  with  thermocouples,  which  required  only  low 
time  resolution,  but  high  sensitivity* 

4.4-2  High  Speed  Photography 

Arc  stability  problems  were  investigated  with  the  Beckman  &  Whitley 
high  speed  cameras  models  326  (Dynafax)  and  307.  The  Dynafax  is  a 
framing  camera  with  a  maximum  framing  rate  of  26,000  pps  and  exposure 
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times  down  to  about  1  microsecond.  The  307  is  a  drum  camera  producing 
"smear"  pictures  with  a  time  resolution  as  short  as  100  nanoseconds.  For 
this  project  both  cameras  were  used  only  at  intermediate  speeds ,  often 
simultaneously  with  good  success. 

4.4.3  Time  Resolved  Spectrograph 

Our  basic  spectrograph  was  the  Hilger  medium  glass  instrument. 

The  dispersion  of  this  instrument  in  the  near  infrared  spectrum  is  only 
about  100  A  /mm  but  this  is  sufficient  to  resolve  the  strongly  broadened 
infrared  lines  of  oxygen  and  nitroger^and  to  observe  the  effect  of  self¬ 
absorption  on  the  oxygen  7772-7775  A  triplet.  Time  resolution  for  our 
spectral  measurements  was  achieved  by  the  use  of  a  mirror,  Fig.  14, 
rotating  at  135  revolutions /sec.  Timing  of  the  spectroscopic  picture 
within  the  duration  of  the  power  current  pulse  was  controlled  by  an  elec¬ 
tronic  delay  and  was  usually  chosen  to  be  near  the  end  of  the  pulse.  In 
the  optical  system  between  arc  and  spectrograph  two  Kodak  lenses  designed 
for  aerial  cameras  were  used.  Because  of  the  limited  aperture  of  the 
spectrograph  (f:12),  they  were  stopped  down  to  a  fraction  of  their  maximum 
aperture.  The  optical  system  clearly  resolved  a  1/Z  mil  tungsten  ignition 
wire  on  the  slit  of  the  spectrograph. 

Spectra  were  recorded  on  Kodak  IN  plates,  which  are  sensitive  up  to 
9000  A,  and  developed  for  5  minutes  in  1:1  Kodak  Dektol  developer.  The 
unusually  long  developing  time  resulted  in  high  contrast  plates,  but  it  enabled 
us  to  photograph  the  spectrum  of  the  carbon  standard  with  the  same  exposure 
time  as  the  air  arc  by  using  a  larger  spectrograph  slit  width.  The  high' 
pressure  air  arcs  were  photographed  with  a  slit  width  of  .  025  mm  or  .  050  mm 
and  the  carbon-standard  needed  a  .  600  mm  wide  slit. 

Intensity  steps  for  the  absolute  measurements  of  arc  brilliancy  were 
generated  by  geometrical  apertures  placed  in  the  spectrograph  between  the 
first  lens  and  the  prism  in  the  parallel  light  beam.  They  were  made  rect¬ 
angular  with  a  variable  height,  h,  and  a  constant  width,  b.  Fig,  25.  Because 
it  leaves  the  base  of  the  prism  unchanged,  this  procedure  does  not  impair 
spectral  resolution. 

The  carbon  crater  of  McPherson  (Me,  1)  was  used  as  an  intensity  standard. 
There  is  some  conflict  between  Euler's  (Eu.  1)  data  from  this  standard  and 
that  of  Null  and  Lozier  (Nu.  1),  which  are  the  two  best  known  sources.  We 
have  plotted  both  sets  of  data  in  Fig.  26.  Euler's  curve  was  calculated 
using  his  crater  temperature  of  3995  K,  his  absorption  coefficients  as  listed 
in  Table  4.  1,  and  the  black  body  brilliancies  of  Pivovonski  and  Nagel  (Pi,  1), 
Null  and  Lozier's  curve  was  obtained  using  their  value  of  3800  K  for  crater 
temperature  and  a  constant  absorption  coefficient  of  0,98#  Null  and  Lozier's 
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TABLE  4.1 

WAVELENGTH  DEPENDENCY  OF  CARBON  CRATER  INTENSITY 


\  (Microns) 


2  o 

I(erg/cm  -A-ster) 


(Euler) 


(Null  &  Lozier) 


2500  3000  4000  5000  6000  7000  8000  9000 


3455555  5 

4.81  2.  16  1.06  2.  13  2.  88  3.  18  3.  13  2.91 


3445555  5 

3.  19  1.  60  8.  88  1.  93  2.  74  3.  14  3.  17  3.00 


10, 000 


a(Eucler)  0.  706  0.  718  0.736  0.  748  0.756  0.  763  0.  767  0.  772  0.  776 


2.  60' 
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values  are  also  included  in  Table  4.  1*  The  brilliancy  values  disagree  con¬ 
siderably  in  the  ultra-violet  and  in  the  blue  region  of  the  spectrum  but  agree 
very  closely  in  the  vicinity  of  S330  R  where  we  have  measured  continuum 
intensities. 

Carbon  crater  intensity  was  recorded  by  removing  the  high  pressure  arc 
assembly  and  locating  the  carbon  crater  in  exactly  the  same  position.  Except 
for  the  slit  width  of  the  spectrograph  and  the  aperture  before  the  prism,  the 
optical  system  was  unchanged.  However,  the  crater  is  viewed  directly, 
while  the  air  arc  is  seen  through  a  Plexiglas  window  in  the  high  pressure 
vessel  so  the  transmission  of  this  window  {calculated  to  be  90%)  was  taken 
into  account  in  our  ^evaluations.  The  crater  was  operated  closely  to  its  hissing 
current  and  in  a  90  mounting, 

4.4.4  Total  Radiation  Calorimeter 

The  calorimeter  developed  for  total  radiant  heat  flux  measurement  is 
shown  in  Figs.  27  and  28.  The  unit  is  a  temperature  sensing  element 
surrounded  by  insulation,  except  for  the  surface  exposed  to  the  arc,  and  sand¬ 
wiched  between  two  copper  discs  of  the  type  used  to  make  up  the  remainder 
of  the  arc  cage  so  that  the  calorimeter  could  be  installed  as  an  integral  part 
of  the  arc  cage.  The  sensing  element  is  a  copper  ring,  to  which  is  welded 
a  copper-constantan  thermocouple. 

As  indicated  in  Fig.  28,  three  designs  were  built.  The  first  two  were 
to  determine  the  relative  effect  of  dishing  the  copper  disc  covers*  The  third 
represented  an  attempt  to  increase  surface  absorption  and  reduce  conduction 
but  was  not  tested. 


The  very  large  voltages  across  the  arc  cage  (up  to  1  5  kv  at  ignition) 
precluded  measurements  during  the  shot.  Two  vacuum  switches,  triggered 
by  the  rotating  mirror  which  controls  arc  firing  and  spectrograph  timing, 
isolated  both  leads  from  the  recording  oscilloscope  until  approximately  15  ms 
after  the  arc  was  extinguished,  so  measured  temperature  rise  was  the  result 
of  an  integrated  heat  flux  during  and  after  the  arc  pulse  (700-800  \x  s). 

The  radiant  flux  per  unit  length  may  be  determined  from  measurements 
of  temperature  rise,  A  T  and  time  of  arc  pulse  A  t. 


(q/z) 


rad 


m 


2h 


A  T 
At 


X 


{4.  1) 


The  mass,  m  ,  specific  heat  c  and  half-width,  h,  of  the  calorimeter 
c  e 

ring  are  fixed  for  a  particular  calorimeter  and  £  is  the  ratio  of  radiation 
per  unit  arc  length  actually  reaching  the  calorimeter  to  that  leaving  the  arc. 
The  method  used  to  calculate  this  ratio  is  discussed  in  Appendix  D. 
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Absorptivity,  <X  ,  of  the  calorimeter  sensing  surface  was  made  close  to 
unity  by  blackening  and,  on  model  3,  threading  to  get  the  effect  of  a  wedge 
shaped  cavity* 

Several  other  corrections  to  the  thermocouple  reading  have  been  considered. 
Applying  the  conduction  equation,  Eq.  (2*41),  from  r  to  the  calorimeter 
radius  yields  a  thermal  conduction  heat  flux  per  unit  length  from  the  arc 
of  the  order  of  1  kw/cm,  but  this  must  be  considerably  reduced  by  the  relatively 
cool  copper  disc  covers,  especially  with  the  dished  configuration  Models  1 
and  3.  Energy  in  the  arc  heated  air  after  power  is  cut  off  was  estimated 
from  enthalpy  and  density  and  is  only  a  few  percent  of  the  total  energy  dissi¬ 
pated* 

An  estimate  of  response  time  was  made  using  a  simplified  model  of  the 
calorimeter  with  the  following  assumptions: 

1) .  During  the  arc  pulse  the  radiant  flux  is  uniform  and  constant  and 

the  calorimeter  sensing  element  is  a  semi-infinite  solid. 

2) p  After  the  arc  is  extinguished  the  sensor  is  completely  insulated  and 

the  heat  flux  within  is  1 -dimensional* 

The  temperature  distribution  at  time  of  arc  extinction,  t  ,  was  computed 
using  the  following  equation  based  on  a  solution  given  by  Carsiaw  and 
Jaeger  (Ca.  !)• 


T 


(x, 


_  2  (q/A) 
1}  k 


<-x2/4a  ctj) 
e 


(4.  2) 


Temperature -time  history  was  then  computed  with  a  standard  finite 
difference  method,  again  making  use  of  the  small  digital  computer.  Results 
are  shown  in  Fig.  29.  The  semi-infinite  solid  approximation  is  clearly 
valid  during  arcing  since  temperature  rise  at  the  outside  radius  is  negligible 
at  time  t  .  For  a  radiant  flux  density  of  5  kw/cm  at  the  calorimeter  sur¬ 
face  (a  representative  value)  thermal  equilibrium  is  approached  to  within 
0.  1  K  in  slightly  over  5  milliseconds.  The  thermocouple  was  located  near 
the  mean  radius  of  the  calorimeter  and  in  that  position  apparent  equilibration 
is  even  faster.  This  response  time  is  substantially  shorter  than  the  time 
needed  to  switch  in  the  thermocouple  circuit. 


The  following  estimates  were  made  of  heat  loss  through  the  surfaces 
in  contact  with  the  calorimeter  ring.  Consider  first  the  Plexiglas  disc  which 
surrounds  the  outer  surface.  Assume  that  the  temperature  of  the  copper 
ring  is  constant  and  that  the  Plexiglas  is  a  semi-infinite  solid  in  perfect 
contact.  Then,  at  any  time  t,  the  total  amount  of  heat  transferred  through 
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unit  area  from  time  zero  can  be  found  by  a  spatial  integration  of  the  temperature 
rise  in  the  semi-infinite  solid.  The  latter  is  well  known  (e.  g.  Ref.  Ro.  1). 

This  heat  must  be  supplied  from  the  copper,  and  the  approximate  fractional 
change  in  calorimeter  temperature  is 


<AT)C 

TC'TP,» 
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2  2 
r  -  r, 
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(4.  3) 


Where  subscripts  C  and  P  refer  to  calorimeter  ring  and  Plexiglas. 
With  property  values  from  Table  4.  2  the  fractional  change  in  calorimeter 
temperature  is 


(AT)c 

Tc  -T 


p,® 


0.22  \/  t(sec) 


(4.4) 


which  amounts  to  a  maximum  possible  error  of  2. 7%  at  the  end  of  15  ms. 


The  effect  of  heat  transfer  through  the  flat  surfaces  is  somewhat  less 
certain  since  it  is  dependent  on  the  heating  of  the  copper  discs  which  form 
the  covers  of  the  calorimeter  assembly.  Ignoring  heat  transfer  to  the  walls 
of  the  pressure  vessel,  the  total  energy  dissipated  in  the  arc  (approximately 
300  joules)  would  cause  an  average  temperature  rise  of  about  9°  K  in  the 
cage  and  surrounding  air.  This  is  the  same  order  rise  as  expected  in  the 
calorimeter  ring  so  the  driving  temperature  difference  should  be  small* 
Comparing  the  heat  transfer  through  the  .  02  inch  thickness  of  alumina  to 
that  needed  to  cause  a  temperature  change  (AT)^  in  (by  symmetry)  half  tho 
calorimeter  ring  yields 


2 a  h  PrC- 


C  C 


(4.5) 


where  z  is  the  thickness  of  the  aluminum  layer  and  h  is  again  the  half  thick¬ 
ness  of  the  calorimeter  ring 


<AT)c 

Tc  -T 


4.  0t{ sec) 


(4.6) 


Assuming  that  Tq  -  T  remains  of  the  same  order  as  the  net  temperature 
rise  in  the  calorimeter,  this  amounts  to  an  error  of  approximately  6%  at  the 
end  of  15  ms.  Thus  heat  exchange  with  insulation  has  a  relatively  minor 
effect  on  the  calorimeter  during  the  time  it  takes  to  switch  on  the  thermocouple, 
after  which  direct  observation  can  be  made  of  temperature  change  with  time. 
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Conduction  at  the  exposed  surface  both  during  and  after  arcing  must  be 
accounted  for  as  discussed  previously.  However,  re -radiation  from  the 
exposed  surface  is  completely  negligible  because  of  the  relatively  small  tem¬ 
perature  rise  involved. 

TABLE  4.2 

SELECTED  PROPERTY  VALUES  USED  IN 
CALORIMETER  CALCULATIONS 


Material 

Density 

Thermal 

Specific 

Thermal 

(g/cm  ) 

Conductivity 

Heat 

Di^usivity 

(watt/cm°K) 

Joule  /g°K 

cm  /sec 

Copper^ 

8.  95 

3.  80 

.386 

1.  12 

Plexiglas 

1.2(2» 

.0021*2* 

1.47(3> 

.0012 

Alumina 

.  0017^ 

Air  ( 300°K,  1  atm)(1> 

.  0011 

.  0028 

1.0 

.24 

(1)  Ro.  1 

e  (2)  Ms.  1 

Sources:  (3)  Mn_  j 

(4)  Estimated  from  values  given  in  Wh.  1 
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5*  METHODS  OF  MEASURING  TEMPERATURE 

Because  of  the  high  pressure  and  large  electric  current  the  arc 
was  expected  to  be  optically  thick,  though  not  M black"  so  that  those 
methods  depending  on  an  optically  thin  plasma  would  not  be  applicable. 

For  this  reason  we  contemplated  use  of  Bartel's  method  (Ba,  1)  which 
depends  on  the  phenomenon  of  self  reversal  of  a  spectral  line  in  an 
optically  thick  plasma  with  a  temperature  gradient.  However,  the  arcs 
at  high  pressure  were  smaller  in  diameter  than  expected,  thereby  re¬ 
ducing  their  optical  depth  and,  while  self-reversal  of  the  oxygen  7772- 
7775A  triplet  was  observed  at  400  amperes,  use  of  Bartel's  method 
appeared  to  be  marginal. 

On  the  other  hand  the  optical  depth  of  the  continuum  was  also  re¬ 
duced  to  the  point  where  it  could  be  used  for  temperature  measurement. 

A  succession  of  difficulties  with  arc  stability  delayed  the  spectr ogr aphic 
results  and  an  interim  means  of  estimating  temperature  was  also  employed 
using  electrical  measurements  and  theoretical  values  of  electrical 
conductivity  as  discussed  below. 

5.1  CONTINUUM  RADIATION  INTENSITY 

The  wavelength  selected  for  determination  of  temperature  from  the 
absolute  intensity  of  continuum  radiation  is  in  the  near  infrared  at  8330 
Continuum  radiation  has  frequently  been  used  for  temperature  measure¬ 
ments  in  both  Low  pressure  ( M.  3)  and  high  pressure  arcs  {Bg.  1,  Bo.  1 , 

Mi.  1). 

We  have  chosen  as  a  theoretical  basis  the  calculations  of  Nardone, 
Breene,  et  al  {Na.  1).  These  calculations  include  free-free  and  free- 
bound  radiation  of  the  electron  in  the  field  of  an  ion,  the  free-free 
radiation  of  electron  and  atom  and  the  O  minus  radiation.  However,  the 
amounts  of  free-free  and  of  free-bound  radiation  obtained  by  Nardone, 
Breene,  et  al  were  multiplied  with  the  slightly  temperature  dependent 
Gaunt  factors  determined  by  Karzas  and  Latter  (Ka.  1).  This  was  a 
correction  of  about  20%.  To  compare  these  results  with  those  of  other 
authors,  we  also  calculated  the  intensity  of  the  continuum  radiation 
using  Unsolds  formula  (Ma.  2),  again  with  the  Karzas  and  Latter  gaunt 
factors.  The  result  is  plotted  in  Figure  30,  along  with  the  Nardone 
curve.  The  two  curves  agree  fairly  well,  within  250  K  at  15,000  K 
for  example. 

Besides  the  uncertainty  in  the  theoretical  results,  our  measurement 
of  the  arc  temperature  from  the  continuum  intensity  is  also  affected  by 
self-absorption  in  the  arc  column.  However,  the  optical  depth  of  the  arc 
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at  the  wavelength  selected  for  measurements  is  considerably  smaller 
than  unity  so  self-absorption  can  be  corrected  for  by  a  rather  simple 
method.  First,  all  evaluations  are  done  with  the  assumption  of  an 
optically  thin  plasma.  This  leads  to  a  first  approximation  for  a  radial 
intensity  distribution,  using  a  standard  Abel  inversion  technique  (Cr,  1). 
From  the  resulting  temperature  distribution  and  the  known  values  of  the 
local  brilliancy,  the  absorption  profile  of  the  arc  is  obtained.  This 
profile  can  then  be  numerically  integrated  for  several  values  of  axial 
distance  to  obtain  total  absoprtion  along  each  path.  Radiation  originating 
in  the  arc  will,  on  the  average,  see  approximately  half  the  total 
absorption.  We  have  therefore  corrected  the  originally  measured  I  (x) 
distribution  for  1/2  of  the  integrated  absorption. 

Such  a  simple  correction  of  the  absorption  effect  in  the  continuum 
can  be  made  because:  1)  the  total  absorption  was  small  at  the  frequency 
measured,  2)  the  continuum  radiation  depends  very  sensitively  on  temp¬ 
erature,  Therefore,  the  correction  of  the  absorption  effect  leads  to 
only  a  minor  temperature  correction. 

5.2  ARC  TEMPERATURE  FROM  VOLTAGE  GRADIENT 

A  preliminary  estimate  of  arc  temperature  at  500  amperes  was 
obtained  from  voltage  gradient  measurements  and  observations  or  arc 
diameter.  An  enthalpy  profile  was  assumed  of  the  form 


where  subscript  a  here  refers  to  the  visible  edge  of  the  arc.  Weber's 
correlation  (We.  1)  Eq.  (2.  50),  for  electrical  conductivity  as  a  function 
of  enthalpy  and  pressure  was  then  used  with  Eq.  (5.  1)  to  obtain  electrical 
conductivity  as  a  function  of  radius.  Then,  writing  ohms  law  in  the  form 
of  an  overall  conductance 

_ 1 

V  r  2  E 
a  avg 

the  left  hand  side  was  determined  from  experimental  measurements  of 
current,  voltage  gradient  and  optical  diameter  and  the  equation  solved 
graphically  for  centerline  temperature  by  plotting  the  overall  conductance 
as  a  function  of  centerline  temperature. 
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6.  ARC  STABILITY 

The  achievement  of  arc  stability  at  a  pressure  of  100  atmospheres 
turned  out  to  be  the  crucial  difficulty  of  the  experimental  part  of  this  study. 
This  was  particularly  so  with  the  copper  constrictor,  which,  for  most  of 
this  work,  was  chosen  over  the  ablation  type  constrictor  because  it  per¬ 
mitted  the  investigation  of  air  arcs. 

Our  first  experiments  at  high  pressures,  with  copper  constrictors 
having  a  hole  diameter  of  5  mm,  showed  strong  instabilities  and  a  much  con¬ 
tracted  arc  column,  Fig.  31.  An  impression  of  the  decrease  of  arc  diameter 
with  pressure  may  be  conveyed  by  Fig,  32  in  which  arc  diameters,  as 
measured  from  the  bright  arc  cores  in  high  speed  photographs,  are  plotted 
against  pressure.  Considerably  improved  arc  stability  was  observed  when 
constrictor  diameter  was  reduced  to  2.  5mm  and  2mm  but  a  kink  type  of 
instability  Could  never  completely  be  eliminated. 

We  went  to  some  length  to  prove  that  no  simple  experimental  side  effect 
caused  the  remaining  arc  motions: 

1)  Current  flow  in  the  arc  apparatus  was  completely  symmetrized  by 
using  the  high  pressure  vessel  as  a  back  conductor.  This  measure  minimized 
undesired  magnetic  forces  on  the  arc. 

2)  The  arc  column  was  made  rather  long,  7  cm,  to  minimize  electrode 
effects  on  the  remote  parts  in  the  middle  between  the  electrodes. 

3)  The  copper  discs  were  judiciously  cleaned  and  etched,  as  described 
earlier,  so  that  only  a  solid  copper  wall  remained  without  loose  particles  of 
copper  or  crystals  on  the  surface.  Such  particles  w'ould,  under  the  heat  load 
of  the  arc,  evaporate  at  considerable  shorter  times  than  the  solid  copper  wall 
and  blow'  at  the  arc  in  an  uncontrollable  manner. 

4)  The  1/2  mil  tungsten  ignition  wire  was  strung  straight  and  in  the 
very  center  of  the  constricting  channel  because  it  was  observed  that  the  arc 
tends  to  remain  very  close  to  the  shape,  in  which  it  was  when  the  wire 
exploded  for  a  rather  long  time  {several  hundred  microseconds)  after  ignition. 
When  started  straight  and  in  the  very  middle  of  the  constrictor,  the  arc 
remains  in  that  position  for  some  time. 

All  these  measures,  however,  did  not  prevent  the  arcs  in  the  copper 
constrictor  from  eventually  going  unstable  in  a  kink-like  form. 

As  for  the  degree  of  instability,  there  was  a  difference  between  the 
400  amp  and  the  100  amp  arcs  burned  at  100  atmospheres  in  2.5  mm  con¬ 
strictors.  The  400  amp  arcs  were  never  really  stable  and  shov’ed  motions  at 
least  in  their  outside  parts.  This  is  demonstrated  with  simultaneous  high 
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speed  photographs  from  a  framing  camera  and  a  smear  type  camera  in  Fig. 
33.  At  lower  currents,  the  arcs  can  be  operated  for  longer  periods  of  time 
without  melting  the  copper  discs.  Then  three  stages  of  arc  development  can 
be  distinguished,  as  seen  in  detail  in  Figures  34,  35  and  36.  During  stage  A 
the  wire  explodes  and  the  arc  plasma  expands  until  it  fills  the  constrictor. 
Then  the  arc  remains  rather  stable  in  stage  B  as  a  diffuse  plasma  column 
filling  the  channel.  Finally  in  stage  C,  a  distinct  thin  arc  core  appears, 
which  moves  slowly  in  the  constrictor.  As  seen  from  the  discrete  pictures 
of  Fig.  36  complete  stability  is  never  achieved  with  the  deviations  from  a 
cylindrical  distribution  growing  in  time. 

We  cannot  say  at  this  time  what  the  reason  is  for  the  very  charac¬ 
teristic  formation  of  an  arc  core.  However,  it  appears  that  the  diffuse 
period  of  the  arc  is  caused  by  the  contamination  of  the  plasma  with  the  vapor 
of  the  1/2  mill  tungsten  ignition  wire.  This  problem  requires  further  work. 

In  addition  to  the  experiments  using  copper  constrictors,  a  few  tests 
were  made  with  ablation  type  arc  channels.  The  channels  were  made  of 
Plexiglas  (Fig.  16}  with  an  inside  diameter  of  3  mm.  The  current  was 
400  amperes  and  the  pressure  150  atmospheres.  Fig.  37  shows  these  arcs 
to  be  remarkably  stable  after  only  a  brief  period  of  instability  immediately 
following  arc  ignition. 

The  difference  between  arcs  in  the  copper  constrictor  and  those  in  the 
Plexiglas  constrictor  is  quite  remarkable.  All  arcs  in  the  copper  con¬ 
strictors  eventually  go  unstable,  even  if  they  were  stable  initially.  Arcs  in 
Plexiglas  constrictors  show  just  the  opposite  behavior  and  become  stable 
after  a  short  initial  unstable  time.  It  appears,  therefore,  that  in  the  copper 
constrictor  the  destabilizing  forces  prevail  and  in  the  Plexiglas  constrictor 
the  stabilizing  forces  are  stronger. 
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7.  RESULTS  AND  CONCLUSIONS 

As  described  in  the  last  chapter*  only  the  ablation  type  arcs  were  com¬ 
pletely  stable*  However,  most  of  our  measurements  were  made  on  arcs  in 
copper  constrictors  which  were  not  completely  stable* 

7.1  ELECTRICAL  GRADIENTS  AND  PRELIMINARY  ARC  TEMPERATURE  DETERMINATION 


Electrical  gradients  were  determined  from  measurements  of  the  total 
arc  voltage  and  the  length  of  the  constrictor.  The  sum  of  cathode  and  anode 
fall,  was,  in  all  cases,  estimated  to  be  50  volts.  At  1500  psi  the  assumption 
of  a  linear  potential  distribution  along  the  arc  column  was  checked  by  using  4 
of  the  copper  discs  of  the  cascade  as  probes.  This  measurement  confirmed 
the  linear  potential  distribution. 

Voltage  gradients  of  500  ampere  air  arcs,  burning  in  copper  constrictors 
of  2,5  mm  diameter  and  7  cm  length,  are  plotted  as  a  function  of  pressure  in 
Fig,  32.  Using  these  gradients  and  the  estimates  of  arc  diameter  also  shown 
in  Fig,  32,  Eqs,  (2.48,  5.  1  and  5.2)  were  applied  as  discussed  in  Section  5,2 
with  results  shown  in  Fig.  38.  Three  values  of  profile  parameter  are  shown 
in  Fig,  38;  x-*  ®  corresponds  to  uniform  temperature  and  x  -  2  corresponds 
to  a  parabolic  profile.  The  bars  at  100  atmospheres  indicate  the  effect  of 
varying  Ta  from  8000°K  to  10,  000°K. 

Since  neither  Ta  nor  the  profile  shape  were  known,  precise  values  of 
centerline  temperature  could  not  be  determined,  but  the  following  tentative 
conclusions  were  reached: 

1)  Since  a  uniform  temperature  profile  is  certainly  a  limiting  case, 
centerline  temperature  is  at  least  16,000  K  for  a  500  ampere  air  arc  at 
100  atmospheres  and  is  more  likely  between  18,000  K  and  20,000  K, 

2)  Temperature  changes  only  slowly  with  pressure,  and  appears  to 
have  a  minimum  value  somewhere  between  20  and  100  atmospheres. 

7.2  TEMPERATURE  MEASUREMENTS  FROM  CONTINUUM  INTENSITIES 

Measurement  of  the  radial  temperature  distribution  were  made  on  arcs 
with  a  pressure  of  103  atmospheres  and  currents  of  115  and  400  amperes, 

7.2.1  The  115  Ampere  Arc 

The  result  of  the  measurements  at  the  115  ampere  arc  is  shown  in  Fig. 
39-  An  axis  temperature  of  17,650  K  was  obtained  and  a  radial  distribution 
with  "core  formation,  11 
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This  measurement  w as  made  close  to  the  end  of  a  2300  microsecond 
current  pulse  where  a  rather  stable  arc  core  had  formed,  as  demonstrated 
by  Fig.  40.  This  figure  shows  a  "smear11  picture  of  the  arc  core  during  the 
spectroscopic  exposure  time.  The  photograph  shows  that  the  core  did  not 
burn  in  the  very  center  of  the  arc  constrictor  and  this  fact  was  confirmed  by 
a  9light  asymmetry  in  the  photometer  curve  of  the  radial  distribution  of  the 
brilliancy.  The  curve  was  folded  and  smoothed  for  the  further  evaluation. 

Fig.  40  also  shows  the  infrared  spectrum  of  the  arc  with  strong  and  very 
broad  oxygen  and  nitrogen  lines  and,  in  addition  to  these,  a  continuum  of  con¬ 
siderable  strength.  QThe  temperature  measurement  was  made  at  the  continuum 
wavelength  of  8330  A. 

7.2,2  The  400  Ampere  Arc 

For  the  400  ampere  arc,  no  reasonable  radial  temperature  distributions 
were  obtained.  This  arc,  as  shown  in  Fig.  33,  fluctuates  strongly  in  the  con¬ 
strictor,  which  tends  to  lower  the  spectroscopically  measured  axis  temperature. 
Consistent  with  this  argument,  an  axis  temperature  of  only  16,  600  K  was 
measured. 

It  was  also  observed  that  the  radial  distribution  of  the  continuous  radiation 
extended  well  beyond  the  wall  of  the  copper  cage.  There  are  several  possible 
explanations,  including  a  penetration  of  arc  plasma  into  the  area  between  the 
copper  discs.  Further  experimental  work  is  needed  to  clarify  this  point. 

7.3  APPLICATION  OF  ARCRAD  lit  TO  MEASURED  PROFILE 

Using  the  least  squares  procedure,  Eqs.  (2.38  and  2.39),  a  temperature 
profile  of  the  form  of  Eq.  2.  35  was  fitted  to  the  measured  profile  as  shown  on 
Fig.  41.  The  fitted  profile  has  Ta  =  12,700°K,  x  =  1.  1343  and  a  =  44°K, 
an  average  voltage  gradient  of  140  volts /cm,  <7  =26  volts,  ana  a  current 

of  111  amperes.  These  compare  with  measured  values  of  160  v/cm  and  115 
amperes,  respectively.  The  results  are  the  correct  order  of  magnitude  and 
can  therefore  be  considered  reasonably  successful  but  they  do  illustrate  the 
difficulty  of  computing  without  experimental  guidance. 

As  an  illustration  of  the  radiant  energy  redistribution  predicted  by  the 
ARCRAD  III  analysis,  Fig.  42  show9  cumulative  net  radiation*  computed 
using  the  arc  profile  of  Fig.  41.  The  curves  are  normalized  by  the  total 
radiant  flux  density  at  the  arc  boundary  {8.  36  kw/cm)  and  are  plotted  for 
three  frequency  ranges  as  well  as  for  the  sum  over  all  frequencies.  At  a 
radius  of  .05  cm,  vacuum  u.  v.  radiation  is  over  25%  of  that  which  eventually 
leaves  the  arc,  but  practically  all  of  the  vacuum  u.v.  is  reabsorbed  in  the 
outer  layers  so  that  effective  emission  is  all  below  59,000  cm’^. 


*  By  Cumulative  net,f  at  any  radius,  r,  is  meant  the  net  radiant  flux 
crossing  a  cylindrical  surface  of  radius  r. 
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7.4  TOTAL  RADIANT  HEAT  FLUX  MEASUREMENTS 

All  total  heat  flux  measurements  were  made  at  the  400  ampere  current 
level  where  no  useful  temperature  profile  data  were  obtained.  Thus  the  cor¬ 
relation  of  ARCRAD  III  computations  and  total  heat  flux  measurements  cannot 
be  direct.  Evidence  of  instability  obtained  after  the  calorimeter  measure¬ 
ments  were  made  also  casts  some  doubt  on  the  assumption  of  steady  state 
upon  which  the  analysis  is  based. 

Table  7.  1  summarizes  the  results  of  the  heat  flux  measurements.  In 
view  of  the  uncertainties  regarding  the  arc  itself  no  corrections  have  been 
made  other  than  the  geometrical  one  discussed  in  Appendix  D, 

The  ratio  C  of  radiation  reaching  the  calorimeter  to  radiation  at  the 
arc  boundary  was  computed  using  a  profile  estimated  to  be  representative. 
This  ratio  appears  to  be  quite  insensitive  to  profile  parameter  as  can  be  seen 
from  Table  7,  2  which  shows  results  of  a  set  of  calculations  for  x  =  50.  The 
The  latter  value  was  chosen  to  reduce  the  number  of  annuli  in  the  calculation 
since  it  could  only  be  a  rough  approximation  in  any  case  and  also  to  render 
the  choice  of  Ta  unimportant.  The  radiant  fluxes  per  unit  arc  column  calcu¬ 
lated  are  of  the  same  order  as  those  measured  although  the  voltage  gradients 
are  somewhat  lower  than  measured. 

Voltage,  current  and  calorimeter  thermocouple  traces  are  shown  in 
Fig.  43  for  the  case  marked  *  in  Table  7.  1.  The  temperature  continues  to 
rise  gradually  as  the  entire  high  pressure  container  tends  to  thermal  equili¬ 
brium.  Traces  taken  with  a  slower  sweep  indicate  that  a  peak  is  reached  at 
about  200  milliseconds  and  after  that  temperature  drops  slowly.  Temperature 
was  measured  as  the  difference  between  a  reference  sweep  prior  to  the  shot 
and  the  temperature  trace  projected  back  to  the  time  of  arc  cutoff. 
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TABLE  7-1 

TOTAL  RADIATION  HEAT  FLUX  MEASUREMENTS  (q/x)rod 


Calori¬ 

meter 

Model 

£ 

Pressure 

(atm) 

Current 

(amp) 

Voltage 
Gradient 
(volts  /  cm) 

<qMrld 

(lew  / cm) 

El 

(kw  /cm) 

hfl 

i 

.24 

20 

115 

9.  5 

49.  5 

.  19 

34 

I 

115 

17.  1 

47.  0 

.  36 

*100 

196 

33.  6 

80.  3 

l  *42 

2 

.  44 

20 

440 

122 

15.4 

53.  6 

.29 

100 

390 

203 

32.2 

79.  0 

.41 

i 

100 

400 

183 

33.  3 

!  73.  0 

.46 

o  o 

Note:  based  on  T^. ^  =  18000  K  Ta=  12500  K  r^  =  .  086  cm 

p  =  100  atm  x  -  1.2  and  2.0  a 


*Fig.  43 


TABLE  7.2 

ARCRAD  III  RADIANT  FLUX  CALCULATIONS  FOR  CONDITIONS 
REPRESENTATIVE  OF  THE  100  ATM  400  AMP  ARC 


Temperature 

°K 

fa 

(cm) 

E 

(volts /cm) 

I 

(amp) 

i  /  \rad 

(q/z) 

(kw/cm) 

c 

(model  1) 

1 

282 

25.2 

.27  ' 

16600 

| 

421 

37.  8 

.27 

.  12 

117 

584 

52.2 

.27 

.  08 

128 

318 

28.  3 

.29 

17000 

.  10 

122 

473 

42.2 

.29 

.  12 

118 

656 

58.  1 

.  30 

.08 

138 

423 

39.0 

.29 

18000 

.  10 

132 

628 

57.8 

.29 

.  12 

126 

867 

79.  0 

.29 
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7.5  CONCLUSIONS 

A  tractable  model  of  the  frequency  dependent  radiant  interchange 
within  a  steady,  cylindrically  symmetric  arc  has  been  developed0  Results 
of  computations  with  the  model  show  that  radiation  in  the  vacuum  ultra¬ 
violet  strongly  influences  the  distribution  of  energy  within  the  arc  and  the 
gas  immediately  surrounding  it.  However,  this  radiation  does  not  contribute 
to  the  wall  cooling  load  in  an  arc  chamber  because  of  very  strong  molecular 
absorption. 

The  capability  for  prediction  of  temperature  profiles  without  some 
guidance  from  experimental  observations  is  limited  because  of  uncertainties 
in  input  radiation  data  {such  as  the  effect  of  lines),  and  in  thermal  and  elec¬ 
trical  conductivities  at  high  pressures.  These  uncertainties  led  to  retention 
of  a  single  parameter  temperature  profile  which  matched  the  measured  profile 
quite  well  and  gave  local  voltage  gradients  of  the  correct  order  of  magnitude 
but  did  not  justify  a  more  detailed  solution. 

The  relative  significance  of  line  radiation  from  a  high  pressure  arc 
might  be  evaluated  by  examining  a  few  typical  lines  in  some  detail  (i,  e. 
estimating  or  measuring  line  broadening  effects  and  taking  perhaps  10  wave 
number  intervals  over  several  half-breadths).  Direct  experimental  verifica  - 
tion  would  also  be  possible  by  measuring  average  intensity  on  a  plate  on  which 
the  lines  are  deliberately  smeared  by  using  a  large  slit  and  comparing  the 
results  to  measurements  of  continuum  intensity  only  from  a  sharp  plate. 

Constricted  air  arcs  have  been  studied  under  conditions  where  very 
little  quantitative  experimental  work  has  been  done.  Arcs  with  currents  from 
100  to  500  amperes  have  been  run  at  100  atmospheres  pressure.  Wall  stabiliza¬ 
tion  was  not  effective  for  maintaining  a  stable  arc  at  such  high  pressure. 
However,  at  the  lower  current  level  of  100  amperes,  a  meta-stable  state  was 
maintained  in  the  air  arc  for  sufficient  time  to  expose  a  spectrographic  plate 
from  which  a  temperature  profile  was  determined.  Such  a  plate  also  contains 
a  great  deal  of  additional  information  on  spectral  distribution  of  continuum 
and  line  radiation  intensities.  Much  additional  data  reduction,  including  use 
of  an  Abel  inversion  which  accounts  for  self -absorption  (Gr.  1,  Bo.l)  will  be 
needed  to  extract  this  information.  Radiant  heat  flux  measurements  are  in 
reasonable  agreement  with  calculations  in  view  of  the  uncertainty  in  arc  state 
at  the  100  atm,  400  amperes  condition  where  heat  flux  measurements  were 
made. 

Arcs  which  are  flow  stabilized  by  abLation  from  the  constrictor  wall 
are  stable  even  at  150  atmospheres  and  prpbably  higher.  Although  the 
resultant  plasma  is  a  mixture  corresponding  to  the  wall  composition  (typically 
some  combination  of  C,  H  and  O)  rather^han  air,  the  inherent  stability  of  the 
ablating  constrictor  type  arc  offers  the  proepect  of  studying  both  local  spectral 
radiant  intensities  and  total  radiant  flux  density  over  a  broad  range  of  pressure 
and  arc  current.  It  is  recommended  that  future  work  on  high  pressure  plasma 
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radiation  place  primary  emphasis  on  the  ablating  constrictor  type  arc 
because  of  its  much  greater  stability  and  flexibility  of  operation. 
Theoretical  data  on  continuum  radiation  of  the  atomic  species  are  available 
and  the  importance  of  line  radiation  could  be  assessed  as  described  above. 
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APPENDIX  A 

TABULATION  AND  INTERPOLATION  OF  TRANSMISSION  INTEGRAL 


An  integral  arises  in  the  radiant  transfer  model  used  in  this  work 
which  cannot  be  solved  in  closed  form,  namely 


f(x)  = 


-W  2 

■f  ■ 


-x/sin(5  .  2  - 

e  sin  0d0 


(A.  1) 


The  original  computer  program  ARCRAD  (Ma.  1,  Ma,2)  took  advantage  of  the 
close  resemblance  of  the  integral  to  an  exponential  decay  by  using  an  empirical 
representation  which  was  the  product  of  an  exponential  decay  and  a  polynominal. 
Since  Program  ARCRAD  III,  which  includes  angular  variation  of  emission, 
requires  the  difference  between  two  values  of  the  integral,  a  more  precise 
evaluation  was  needed.  The  function  was  integrated  numerically  to  an  error 
less  than  1,0  x  10"^  using  the  parabolic  (Simpson's)  rule.  A  table  of 
differences  was  then  calculated  and  stored  in  the  computer  to  facilitate 
evaluation  of  the  integral  by  means  of  Bessel's  interpolation  formula  (Hi.  1) 


f(x)  =  d  (x.)  +  (s-1/2)  d  (x.)  +  s(s-l)  d  (x.) 

O  1  1  1  Z  1 

+  s(s-l)<s-l/2)  d3(x.) 


where  x.  is  the  next  lower  value  for  which  x  is  tabulated  and 
i 


x  -  x. 


s  = 


x  -  X. 

1+1  1 


«*.)  +  f(*i+1> 
d  (x.)  =  - = - 


O  1 


d^Xf)  =  f(x.  +  1)-f(x.) 


dJx.)  = 

Z  i 


£t*l+2>  -f(xi+l>  -«*i>  +£(’ti.1) 
2  x2! 


d3(x.)  = 


£(X.+2)  -  3f(*.+1)  +  3f(x.)  -  £(Jtt  l) 
2x3! 


(A.  2) 

(A.  3) 

(A.  4) 
(A.  5) 

(A.  6) 

(A.  7) 
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The  interval  was  varied  with  optical  depth  to  keep  the  table  of  differences 
as  small  as  possible  while  retaining  a  precision  of  interpolation  equaling 
the  precision  of  the  tabulated  values  of  the  function* 

Table  A.  1  lists  tabulated  values  of  f(x)  and  associated  differences  for 
x  =  0  to  x  =  14*  At  larger  optical  depths  interpolation  of  ln(f(x))  proved  more 
convenient,  so,  in  the  range  6  ^  x  ^  14,  the  table  lists  In  (f(x))  and  its 
associated  differences.  At  an  optical  depth  of  14,  f(x)  <  .5  x  10"°  so  this 
was  chosen  as  the  cutoff.  Fortran  compatible  notation  was  used  in  the  table 
with  F(NX),  dQ(xi)^F(NX,  1),  etc. 
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table  A*1 

TRANSMISSION  INTEGRAL  AND  INTERPOLATION  DIFFERENCES 


ft  KX )  F(NX4)  F(NX,2}  FJNX,3] 


hX 

■(NX) 

*TNXfl  1 

F(NXt  2) 

FI  NX,  1J 

F(NX.  4 1 

.20 

+  7. 793^3^-*  t 

+7,5B20BBv?E-01 

-4.6191  46PE -02 

>1 .56528E-C3 

.00 

>9.9999946f-WI 

+9.87461 7tfE-01 

-2.5P7543PE-G2 

+3.66450E-04 

-3.791 86E-06 

.26 

+  7.33M  290£-f  1 

+7.M51690E-01 

-4.3191 9  70E-02 

>1,44178E-03 

.02 

>9.749z403E-0 i 

♦9.6276135E-01 

-2,432533«f-fl2 

+3.66672E-4I4 

-3.13833E-06 

.60 

+6.B992093E-«1 

+6.697PB75E-01 

.0424330E-02 

>1 . 331 97E — 03 

.04 

♦  9.5©5937<3E-fl  1 

♦9.3B79230F-01 

-2.361274HE-02 

>3.43232E-0'« 

-2.68SB3E— 06 

.35 

>6.4949660E-01 

+6.3056455E-0M 

-3.7B640  70E-«2 

+1.23346e-*3 

. at, 

>9.269B59&E-*I 

>9.155I975E-01 

-2'.29324«lE-02 

>3.33PI0E-04 

-2.*39I66E-«6 

,40 

+£.11 63253EH91 

>5.93687  25E-A l 

-3.5490470E-02 

>1 .1 444> E— 03 

.OR 

>9.04P5  356E-«!1 

>fi,929132PE-0l 

-2.228O7ME-02 

+3  .*19375E-04 

-2.  16333E-4J6 

.  43 

+5.7B14206E-01 

>5  .5949890E-01 

-3.32B6260E-02 

+  1  .063S5E-03 

.1  o 

77286 F-fM 

>8,7094540E-41 

-2',16549m-«2 

+  3.^696CE-C4 

-1.9BP00E-0G 

.50 

+5>2e5?80E-d1 

+5.27  237Fi5E— C 1 

-1,1 236270E-0  2 

4-9.B9782F-24 

At 

tB.fi{in796E-f)l 

>8.49591 50 E-01 

-2.1d52860E-{'2 

+  2.954$7E-flf* 

-1  .855tWF-06 

+5 .1 1^79^3  t-fl 1 

iii.9h9559^E-P  1 

-2.93271  3C1E-02 

+9.22257E-04 

.14 

>8,2882a6(iE-fll 

-2*flr*72971E-02 

+2.84772E-04 

-T.72166E-C6 

.60 

+F.fi2.'’92»»£!E-t!1  1 

+4.6Br.  187BE-0  1 

-7.7547240E-02 

+8.602 17E-04 

.16 

+9.1 85921 3E-CM 

i 

>8.08635  2i»F~SH 

-U9913770E-#2 

♦2.74742E-C4 

-1.625P0F-06 

.65 

+4.54745  I6F-01 

>4.4 180203E-0 1 

-2.5886260E-C  2 

+B.03I  I7E-04 

.1  R 

+7.9BG7B36E-01 

>7.8699’' 35E-01 

-T.9374903E-02 

>2.fi53lCE-«4 

-1.52750E-06 

.70 

+4.  2885  69CE-01 

+4.166915 IE-01 

-2.4334770E-02 

+7.5P472E-04 

.20 

>7.7930436E-tM 

|>7.69e7BC5E-«n 

-1 .BS52530E-02 

>2.5&367E~*4 

-l.i»S5(WE-*& 

.75 

{*■4,045  241 3  E-0  1 

>3.930B194E-C1 

-Z.2804370E-0  2 

+7.01785E-0U 

.AO 

f3.8|6397^F-T1 

+3.70B7694E-C1 

-2.1527630E-02 

+6.5669ZE-04 

FINX,  4 ) 
2.1fi55i2E-05 

1 .93133E-05 

I  .72916e-«)5 

1.5552SE-05 

- 1  .'41  1 66E-05 

-1.2960BE-05 

- 1 . 1 73O0E-O5 

-1.07833F-P5 

-9.9O25ae-0ft 

-9.13333E-06 

-8.42I66E-06 

-7.6OB33E-06 

-7.2275CE-06 
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(Continued) 


NX 

KNX) 

F(NX.  i) 

F(NX.  1 > 

FjNX,  7) 

F(NX,  i) 

NX 

I(NX) 

1 

F<NX.  11 

F(NX,  2) 

F(NX,  3) 

F(NX,  4J 

,85- 

<♦3.601 12UE-01 

+3. 4 998 3 33 E-01 

-  2. V!  25  7600  E  “02 

+6,l4925£-0^ 

-6.69666E-06 

1.6 

1 

+  1.22  44  203  E-01 

+  • • 1598 162E-0  1 

-  1  t292k38U3E-02 

+  7.40022E-04 

-1.45758E-05 

.SQ 

+3.398S4S3E-01 

+3.30320568-01 

-  1.9067910F-02 

+5 .76130E-04 

-6.23583E-06 

1.9 

+1.0952122E— 01 

+  1, 037571  IE-01 

-1. 152822CE-02 

+6.57667E-04 

t 

-1. 28766 E-05 

.95 

+3.207&6*0E-«1 

+3.11 81006E-0 1 

-1.7953030E-02 

+5.40052E-C4 

-5.79003E-06 

2.0 

♦  9. 7993006E-0Z 
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l.CO 
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+2.91.37966E-01 

U-1.690772SE-02 
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2.1 
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- 1 ,008 1 3E-05 

2.2 
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+4.63399E-04 
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1*0 

+3.02B3352E-01 

+  2.86416I3EH81 

-3.2834  770E-0  2 

+1.96793E-03 
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2.3 

+  7.03 

+f.,*63Bl  i«E-P2 

-7.33516^-03 

+  4 .  1 2823 F -04 

-7.9244iF-f6 

1.1 

+2.6999875 £-01 

i 

♦  2.554241  IE-0  1 

-2,91 492^fi  £-02 

+I.73385E-03 

-3.62975E-05 

2.1 

+6,297  'R93E-D2 

+5.9692  60 E-C 2 

-6.557057:e-C3 

+3.67953E-04 

-7.03275E-06 

1.2 

♦2.4084947E-0I 

+  2.2W979E-01 

-2.58933&0E-02 

+1 .53P2SE-33 

-3.16533E-05 

2.  S 

+5.64 13536F -02 

+5*  348 186PE-B2 

-5,863346t1E-0  3 

+3.281  1 5  E  — fl  4 

-6.24775E-06 

i 

I.l 

+  2*,  1 1*95^  1 1 E-0  1 

+2.0343547E-01 

-2.3028280E-02 

♦1.35 197E-03 

—2.7  68  8  3  E— SI  5’ 

I 

2.6 

+5.C55019PE-0  2 

+4.7927893E-C2 

-5.244594PE-03 

+  2. 92  720  E -04 

i 

-5 .55  1 33E-06* 

1.1 

♦  1 ,9'92*R3E-fll 

♦  1  .Rlfi75irtE-01 

-2.049H,^e-02 

+1.19609E-03 

-2.42733E-05j 

2.7 

+4.5305596E-C2 

+4.2959364E-P2 

-4.6924630E-C3 

♦2.6I251E-04 

-4.93833E-06 

1.5 

♦•.7l4?93  E-01 

+ 1 ,62V69t  E-0ij 

-1  .R24492  'e-02 

+I.05933E-C3 

-2.13 ioBE-05 

2.8 

+4.061 3 133E— 02 

+3.R513339E-P2 

-4.1995R70E-C3 

+7.33255E-04 

-4.39466E-06 

1.5 

+1.5318  ‘'-5t-0i 

+1. 45n?739t-0l 

-l.oZS^IirE-02 

+  9. 3915  2E-04 

-1.87508E-05 

2.9 

+3.6413546E-02 

+3.4533825E-C2 

-3.75944 I7E-P3 

+  2.08336E-04 

-3  .9 1  275E-06 

1.7 

+1.3693P34E-01 

+  1. 29686 18E— D 1 

-1.44983  ICE-02 

+B.33325E-04 

-1. 65250E-05 

3.0 

+3.26541*51-07 

+3.09709R4E-P  | 

-3. 36 624 20 E- 03 

+1.B6144E-04 
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NX 

'(NX) 

F1NX.  1] 

F(NX,  1 1 

F{NX.  11 

F(NX.  41 

NX 

KNX) 

F<NX.  1} 

F<NX,  1) 

FINX.  31 

FfNX,  4) 

3J 

♦2.9207063F-S2 

♦  2.778*1I*3«F-e2 

-3.* l^065PE- C3 

♦1.66368E-04 

-3.10666J-06  , 

4.4 

♦7.21Z745fE-0l 

♦6. 8469475 E-03 

-7.3l6l6tWE-04 

♦3.9<559E-05 

-7.18858E-07 

3.2 

♦  2/6272990E-C2 

♦  2  .49226 1 4E-02 

-2.700760PE-03| 

♦1.40741E-04 

-2.76908E-0e 

4.5 

♦  4#lpi  139^e-P3 

♦6. 15273 1  H'E-03 

-6.568I700E-04 

♦3.54691E-05 

-6.43466F-07 

3.3 

♦  2.3572230E-02 

♦2.2362280E-02 

-2,4  199000E-03 

♦1.33022E-04 

-2.4706lbE-0b 

4.6 

♦'  .P24?^?^C-03 

♦5,529*‘5?5E-03 

-5.6973960E-04 

♦  3.18O97E-05 

-5.76383E-07 

1.4 

♦2.1 1S2330E-02 

♦2.0067990E-02 

-2.16060ME-03 

♦1.18996E-04 

-2.2041*  1E-0  6 

4.7 

♦5.234$ft3«E-f 3 

♦4, 9697 940 E-0 3 

-5.29578P0E-04 

♦2.95323E-C5 

-5 . 162 16E-07 

3.5 

♦  l .098365PE-0  2 

♦  1 .001 I693E-02 

-1.9439 130  E-0  3 

♦  1 .0b<*78E-04 

-1 .96850l“06 

4.8 

♦4.7A50050E-03 

♦4.467I998E-03 

-4.7561040E-04 

♦2.55966E-C5 

-4.62333E-07 

1.6 

♦  1 ,7f139737E-0  2i 

♦  1  ,61603*^-02 

-1.742760C*E-03 

♦9.5  2992F-C5 

-1.75775E-06 

4.9 

♦4.2293946t-03 

♦4.0 157980E-03 

-4.271916CE-04 

♦2.29663E-C5 

-4. 14416E-07 

3.7 

♦  *.*29*9^-02 

♦  i.4ri«;6i  ie-@2 

-1.562716^-0  3 

♦0.53 175E-P5 

-1.5695*^-06 

5.0 

♦3.8022030E-03 

♦  3.61O3305E-03 

-3.8374500E-04 

♦  2.060  94E-05 

-3.71300E-07 

3.8 

*1.373r*253C-02 

♦  1.3^335^-02 

-1.401  490OE-03 

♦7.63977E-05 

- 1.403  75£“0<- 

5.1 

♦3.4 184500E-03 

*3.246081  IE-0 3 

-3.4475370E-04 

♦1.64969E-05 

“3.32833E-0; 

1.4 

1  •  2  3  32  75  *>  E— 0  2 

♦  1. 170M92E-02 

- 1 .2571 290E-03 

♦6.B4247E-05 

“1.25391E-06 

5.2 

♦3.0737043E-03 

♦  2.9188257E-C3 

-3.09757 1PE-04 

♦  1 ,66P3?e-05 

-2.98425E-07 

4. a 

tl.1075630E-02 

♦1 .051 1 732E-02 

-1.127799OE-03 

♦6.12992E-05 

-1.1 2166E-06 

5.3 

♦  7.7439,‘72E-03 

♦2.624776LE-03 

-2.7034  160E-«1* 

♦  1 .49051  E-05 

-2.67525E-07 

4.1 

►9.9478313E-03 

■t-9.<i  i*  106**5  E-03 

-1,01 19327E-03 

♦5.49242E-05 

-1.00  3B0E-06 

5.4 

♦2.4ff\4i*<E_03 

+2.36(3c374E-f3 

-2.50IJ640E-04 

♦I.33B25E-06 

-2.400PBE-07 

1 

*a 

♦8.9366986E-C3 

♦8.4Bl8i*70E-03 

-9.081P230E-04 

♦4.92230E-C5 

-8  .  974ME-07 

5.5 

♦2.2354692E-03 

♦2.12306J^e-03 

-2.248  1 1 3^E-04 

♦1.20168E-05 

-2.15225E-07 

4.3 

♦8.02779ME-03 

♦7.62027S5E-f3 

-0.150407PE-04 

♦4.t*l  215E-05 

-8.031 33E-C7 

5.i 

♦2.0106579E-03 

♦  I.9096234E-03 

-2.0226890E-04 

♦  1 .0791 9E-C5 

- 1 .93091 E-07 
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NX 

f(NX) 

P(NX,  1) 

F(NX,  2} 

F(NX.  3) 

F(NX,4) 

5*7 

+1.8085&90E-03 

+1.7177672E-03 

-1.8164360E-04 

+9.69290E-06 

-1.73250E-07 

5*8 

+1.6269454E-03 

+1.5452967E-03 

-1.6329730E-04 

+8.70675E-06 

-1.55466E-07 

5.9 

+  1.463648  IE-03 

+1.3902398E-03 

-1 .4681660E-04 

+7.82177E-06 

-1.395 25E-07 

6.0 

+1.3168315E-03 

+  1 .2508264E-03 

-1.3201020E-04 

+7.02747E-06 

-1.25241E-07 

NX 

f(NX> 

Ln(F(NX,  1)> 

Ln(F(NX,  2)) 

Ln(F(NX,  3» 

L,n(F(NX,  4)) 

6.0 

+  1 , 3  T  683 1 5  E— 03 

-7.1593245E-00 

-1.05359b3E-00 

+3.02840E-03 

-1.16450E-04 

7,0 

+4.5915480E-04 

-8.2102420E-00 

-! .0482386E-00 

+2.42970E-03 

-8.30 1 6GE-05 

8.0 

+1.6095886E-04 

-9.2S63000E-00 

-1.0438775E-00 

+1.99615E-03 

-6.14916E-05 

9*0 

+5.6^7149ftE-05 

-1.0298366E+01 

-1.0402540E-00 

+1.67087E-03 

-4.69416E-05 

10,0 

+  2.002^7 1 T  E-0S 

-1  . 13370 90E+01 

- 1.037 1940E-00 

+1.42020E-03 

-3.65666E-05 

11*0 

+7»098,'>7f't5E-06 

-1.2372973E+01 

-1.0345730E-00 

+1.22300E-03 

-2.91666E-05 

1 2,0 

+2.5224984E-06 

-1.0323020E-00 

+1.0 6475 E-03 

-2.35833E-05 

13,0 

+8.9847906E-07 

-1.44377T9E+01 

-1.0303140E-00 

+9.35500E-04 

-1.95000E-05 

14.0 

+3.2066249E-07 

-1. 54671 56E+01 

-1.0285600E-00 

+8.28750E-04 

-1.60833E-05 
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APPENDIX  B 
PROGRAM  ARCRAD  III 


B.l  PURPOSE 

The  purpose  of  ARCRAD  III  is  the  calculation  of  radiation  heat  flux 
from  arc  heated  air  based  on  spectral  emission -absorption  characteristics 
of  air  and  a  measured  or  assumed  temperature  profile.  Included  are  com¬ 
putations  of  flux  reaching  a  total  radiation  calorimeter  and  of  voltage  gradient 
based  on  a  local  energy  balance,  and  total  electric  current. 

B.2  GLOSSARY  OF  TERMS 


Fortran  Name 

Definition 

AO,  Al,  A2,  A3 

Constants  for  polynomial  which  multiplies 
transmission  function  for  calorimeter 
calculation. 

AINT 

Intermediate  value  in  calculation  of  R  array. 

ALPHA  (6) 

Angle b  used  in  calculation  of  sector  radii. 

ANGLE 

Angle  sectors  of  quadrant.  (<rf  +  A  <^/2) 

BCDREC  (25) 

Array  in  which  alphameric  identification  of 
tape  is  read. 

BIG  (20) 

Array  for  storing  largest  absorption  length 
for  each  annulus. 

BINT 

Intermediate  value  in  calculation  of  R  array. 

C2T  (20) 

Array  of  C^/T  used  in  calculation  of  v. 

CEM 

Calorimeter  emission-transmission  integral 
for  particular  value  of  ETKDR. 

CONST 

Average  particle  density  of  outermost  (low 
temperature)  annulus. 

CQRAD  (10) 

Calorimeter  integral  summed  over  wave 
number  group. 

COSNA  (6) 

Cosines  of  ALPHA  array. 
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CTEMP  (21) 

CTR 

C  TRANS 

CURRNT 


Array  of  temperatures  for  which  electrical 
conductivity  is  tabulated. 

Calorimeter  transmission  integral  for 
particular  value  of  TKDR. 

Calorimeter  transmission  function, 
emitting  annulus  to  arc  boundary. 

Current. 


DALPHA 

DEV 

Drv 

DH 

DL 

DSINA  (6) 

DSUM 

E  (21) 
EAVG 
EM 


Angle  difference  in  degrees  between  ele¬ 
ments  of  ALPHA  array. 

Weighted  standard  deviation  of  calculated 
voltage  gradient. 

Length  of  subrange  for  numerical 
integration. 

Relative  cumulative  spectral  radiance  at 
upper  wave  number  limit  of  a  particular 
wave  number  increment. 

Relative  cumulative  spectral  radiance  at 
lower  wave  number  limit  of  a  particular 
wave  number  increment. 

Difference  in  the  sines  of  consecutive 
sector  angles. 

Intermediate  value  for  computing  relative 
cumulative  spectral  radiance  when  v  >  1. 

Annulus  voltage  gradient. 

Weighted  average  voltage  gradient. 

Transmission  integral  for  particular  value 
of  ETKDR. 


EKDR 


Dimensionless  emission  length. 


EKDRMX  (20,  76)  Temporary  storage  array. 


ESUM 


Intermediate  value  used  in  calculation  of 
weighted  average  voltage  gradient. 


55 


A  E  DC-T  R  -66*258 

ETA 

ETKDR 

EX  (10) 

EX1  (10) 

EX2  (10) 

F  (88,  4) 

I 

II 

IF  (10) 

II 

IJK 

IN 

IN  1 
INF 


INN 


Black  body  energy  fraction  within  a 
specified  wave  number  interval  (DH-DL). 

Dimensionless  emission,  transmission 
length. 

Temperature  profile  parameter. 

Temperature  profile  parameter,  Zone  I, 
when  Zone  II  is  omitted. 

Temperature  profile  parameter,  Zone  I, 
when  Zone  II  is  included. 

Array  of  interpolation  differences  for 
dimensionless  absorption  length. 

Wave  number  index  or  annulus  subinterval 
index. 

Lowest  wave  number  at  a  given  temperature 
for  which  v  s  .2. 

Array  of  index  limits  for  breaking  total 
radiation  into  wave  number  groups. 

Subinterval  index  for  evaluating  electrical 
conductivity. 

Index  in  calculation  of  relative  cumulative 
spectral  radiance. 

Index  of  outermost  annulus  for 
ARCRAD  III  calculation. 

IN  +  1 

Consider  a  ray  emitted  from  an  interior 
annulus  surface  passing  through  an 
absorbing  annulus.  INF  is  the  index  of  the 
first  boundary  of  the  absorbing  annulus 
encountered.  NF  is  the  index  of  second 
boundary  encountered.  NF  may  equal 
INF-1  or  NT  +  1  -  INF. 

Index  for  bypassing  Zone  II  data  when 
Zone  I  only  is  being  calculated. 
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IPOS 


Index  on  pressure. 


IPST 


Index  for  locating  absorption  coefficient 
on  an  input  record. 


IR 

ITEST 


IW 

IX 

IZONE 


Index  on  output  radiation  values  stored  by 
wave  number  group. 

When  index  on  absorbing  annulus  greater 
than  NT,  ITEST  =  0;  when  index  less  than 
or  equal,  ITEST  =  1. 

Wave  number  index  used  for  finding  II. 

Index  on  profile  parameter. 

IZONE  a  0  when  Zone  II  is  omitted. 

IZONE  =  1  when  Zone  II  is  included. 


J 


Index  for  temperature  interpolation  of 
conductivity. 


KA 

L 


Index  on  annulus  absorbing  from  interior 
surface  of  emitting  annulus. 

Index  for  reading  title  card. 


N 


Index  on  emitting  annulus. 


N 1 

NA 


Number  of  annuli  exterior  to  an.  exterior 
emitting  surface. 

Angle  index. 


NALPHA  Number  of  angular  subintervals. 

NANG  EE  NALPHA  +  1. 


NB  IN  -  1. 

NC  IN 


NCASE 


NDIV 


Number  of  cases;  each  different  combina¬ 
tion  of  temperature  and  pressure 
constitutes  a  case. 

Number  of  subdivisions  for  integrating 
conductivity. 
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NF 

NIF 

NJ 

NJ2 

NF 


See  INF. 

Number  of  wave  number  groups. 

Index  on  absorbing  annuli  for  exterior 
emission. 

Index  used  in  calculating  annulus  and 
sector  radii. 

Index  used  for  summation. 


NQ 

NS 


NS  1 


Index  used  for  annuli  absorbing  from 
exterior  surface  of  emitting  annuluB, 

1)  Index  used  in  calculating  annulus  and 
sector  radii. 

2)  Index  used  to  select  proper  annulus 
and  sector  radius  for  the  emission  length 
of  the  emitting  annulus. 

NS  +  1. 


NT  Number  of  annuli  of  the  profile. 

NT  1  NT  +  1. 

NT  2  2*  NT. 


NT  21  2*  NT  -  NC. 

NT  22  2*  NT  -  NB. 

NW  Number  of  wave  numbers. 


NWBCD 

NWBIN 

NWI 

NWF 


Control  parameter  (=  0  to  read  entire 
BCD  record). 

Control  parameter  {=  0  to  read  entire 
binary  record). 

First  wave  number  of  wave  number 
group. 

Last  wave  number  of  wave  number 
group. 
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NX 

NX  1 

NX  2 

NXX 

P 

PCNT 

PROP  (20,  76) 
PR  TINT  (20) 

QABS  (20,  10) 

QCON  (20) 

QCONA 

QCOR 

QEMIT  (20,  10) 

QNET  (20,  10) 
QRAD  (10) 

QV  (10) 

R  (40,  40,  3) 
RA 

RAD  (20) 


NX  =  NX1  when  Zone  II  omitted. 

NX  =  NX2  when  Zone  II  included. 

Number  of  profile  parameters  when 
Zone  II  omitted. 

Number  of  profile  parameters  when 
Zone  II  included. 

Index  used  to  read  array  F  (88,  4) 
Pressure. 

Interpolation  fraction. 

Linear  absorption  coefficient. 

Electrical  conductivity  weighting  factor 
(SIGINT(N)  /S1GSUM). 

Radiant  flux  per  unit  arc  length  absorbed 
by  annulus  from  all  other  annuli  {summed 
over  a  wave  number  group). 

Net  thermal  conduction  per  unit  arc 
length  from  an  annulus. 

Thermal  conduction  per  unit  arc  length 
at  boundary  of  Zone  I. 

Total  power  per  unit  arc  length. 

Radiant  flux  per  unit  arc  length  emitted 
by  annulus  ( self -absorption  included). 

QEMIT  -  QABS  for  one  annulus. 

Summation  of  QNET  over  all  annuli. 

Net  radiation  per  unit  volume. 

Array  of  annulus  and  sector  radii, 

Radius  at  boundary  of  Zone  I. 

Array  of  radii  ordered  from  outermost 
annulus  to  core. 
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RDALPH 
REC  (700) 

RL 

RR 

RW 

SIGINT  (20) 

SIGMA  (21) 

SIGSUM 

SIGTA 

SINA  (7) 

SQABS  (20) 
SQEMIT  (20) 
SQNET  (20) 
SUM 

T  (2 1) 

TA 

TAA 

TCL 


Array  in  memory  for  reading  in  tape 
records. 

I 

Inner  radius  of  outermost  annulus  divide^ 
by  RA. 

Radius  subinterval  for  conductivity 
calculation. 

5 

Ratio  of  radii  of  Zone  I  and  Zone  II. 

fr°  ^ 

I  a  rdr 
i 


Input  array  of  electrical  conductivity  at 
1000  K  increments  (pressure  dependent). 


a  rdr 


Electrical  conductivity  at  boundary  of 
Zone  I. 

Sines  of  consecutive  angle  sectors  of  the 
quadrant. 

QABS  summed  over  all  wave  numbers. 

QEMIT  summed  over  all  wave  numbers. 

SQEMIT  -  SQABS. 

Intermediate  value  used  in  several 
summations. 

Intermediate  temperature  for  electrical 
conductivity  integral. 

Temperature  at  boundary  of  Zone  I. 

TA/1000  -  used  in  FINANN  subroutine. 

Centerline  temperature. 
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TOLL 

TCL/1000  -  used  in  FINANN  subroutine. 

TEMP  (20) 

Array  of  annulus  temperatures  ordered 
from  outermost  annulus  in  to  core. 

TEMPRW  (20) 

_4 

TT  0*  T 

TL 

Temperature  of  inner  boundary  of  outer¬ 
most  annulus  {2500  K}. 

TKDR 

Dimensionless  transmission  length. 

TR 

Transmission  integral  for  a  particular 
value  of  TKDR. 

TRANS 

Transmission  function  from  emitting 
annulus  through  absorbing  annulus. 

TRANS  1 

Transmission  function  from  emitting 
annuLus  to  absorbing  annulus. 

TT 

Outer  boundary  temperature  of  annuli  in 
Zone  II. 

TWALL 

Temperature  at  outside  boundary  of 

Zone  II  (wall  temperature). 

V 

C2n/T  =  hv  /kfiT 

V2 

v2 

V  4 

v4 

V6 

v6 

VINC 

Radius  increment  for  conductivity 
calculation. 

WBB  (20,  76) 

TEMPRW  (N)  *  ETA 

WBBR 

WBB  (N,  I)  *  RAD  (N) 

WN  (76) 

Wave  number  array. 

WNS  (10) 

Largest  wave  number  of  wave  number 
group. 
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WVNO  (20) 

XT  (88) 


Fortran  Name 


Wave  number  for  maximum  emission 
length  -  associated  with  given  annulus. 

Dimensionless  absorption  length. 

Subroutine  FINANN 

Definition 


CO,  Cl 

DRAD 

F(I) 

FINTG 

I 

J 

K 

L. 

LL 

N 

NT  J 
NTCLJ 

R(N) 

RAD(J) 

RAD  1 


Constants  for  linear  interpolation  or 
extrapolation  of  annulus  inner  radius 
from  latest  two  trails, 

Radius  increment. 


Integral  function. 

Index  on  Simpson  rule  integration. 

Index  used  for  matching  subroutine 
results  to  main  program. 

Index  for  reshuffling  RI(K)  and  TAVG(K) 
in  preparation  for  new  trial.  RI{3)-*RI(2), 
RI(2)-*RI{1),  etc. 

Iteration  index  L  is  incremented  1  during 
each  iteration  but  also  set  back  one  after 
third  and  subsequent  trials. 

Index  used  for  special  case  of  core  and 
innermost  annulus. 

Index  on  temperature  and  radius. 


Index  on  centerline  temperature. 


Radius  where  temperature  is  T(N). 

Outside  radiuB  of  annulus  in 
centimeters. 

Local  radius  for  integration. 
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RCK 


RCOL 

RI(N) 


Annulus  outside  radius  minus  RMIN,  If 
RCK  does  not  exceed  R(N)  then  T(N)  iB 
skipped  and  calculation  proceeds  to  next 
higher  temperature  level. 

Arc  radius  or  Zone  I  radius,  r  . 

a 

Annulus  inside  radius. 


RMIN  Minimum  permitted  value  of  RO(N)  -  R{N), 

set  at  .  025  r  ;  also  see  RCK. 
a 

RO(N)  Current  value  of  RAD(J).  Both  designations 

used  in  order  of  simplify  matching  of  sub¬ 
routine  in  ARCRAD  III. 


SUM  1  ' 
SUM  2  ' 
SUM  3, 


Intermediate  values  in  integration. 


T(N) 


Temperatures  at  which  data  are  tabulated 
(divided  by  1000  K  for  the  subroutine). 


TAVG(L)  Calculated  average  annulus  temperature 

divided  by  1000°K.  Must  match  T(N)  + 
TTOL  before  going  to  next  inner  annulus. 

TCL  Arc  centerline  temperature. 


TDIFA 


Absolute  value  of  TDIF. 


TDIF 


TAVG(L)  -  T(N),  see  TAVG(L) 


TDIFF 

TEMF(J) 


TCL  -  T{1),  used  in  calculation  of  R(N) 
from  profile. 

T(N)*  1000,  see  ROCM. 


TEX 


Exponent  on  temperature  in  integration, 
set  at  8,  see  Appendix  C. 


TRD  Temperature  at  annulus  outside  radius 

calculated  from  temperature  profile. 

TTOL  Allowed  difference  between  T(N)  and 

TAVG(L),  set  at  .002. 
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TXC  ? 
TXC  lj 

TXX 

X 

Y 


g 

T{NTCL)  jUsed  for  special  case  of  core 
and  innermost  annulus. 

T{1),  value  from  main  program. 

Parameter  on  analytic  temperature 
profile. 

Local  temperature  for  temperature 
average  integration. 


Subroutine  LOOK 


Fortran  Name 


Definition 


AINC 
FX 
FY 

NY 

S 
X 

B.3  INPUT 

1.  Tape 

Binary  tape  contains  linear  absorption  coefficients  for  tem¬ 
perature  273°  and  3000°K- 20,  000  K.  This  is  a  tape  with  19  files, 
one  file  per  temperature. 

2.  Cards 

a.  number  of  cases  (fixed  point) 

b.  cards  per  case 


Intervals  for  which  F  array  is  tabulated. 

Transmission  integral  for  value  of  X. 

Calorimeter  transmission  integral  for 
value  of  X, 

Integer  variable  which  indexes  the  nearest 
tabulated  value  of  F(x)  for  XT  less  than  or 
equal  to  X. 

Interpolation  fraction. 

Dimensionless  absorption  length. 
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1)  identification  card 

2)  NAM  1 

TA  zone  I  boundary  temperature 

HA  zone  I  boundary  enthalpy 

RA  zone  I  boundary  radius 

TWAL.L  temperature  at  wall 
KW  ratio  of  radii  of  zone  II  to  zone  I 

TCL  centerline  temperature 

3)  NAM2 

DALPHA  angle  difference  in  degrees  of  elements 
of  ALPHA  array 

IPOS  pressure  of  the  case  index  dependent  on 

pressure  used  in  picking  off  the  proper 
absorption  coefficient 

SIGTA  electrical  conductivity  at  zone  I  boundary 

4)  NAM3 

NIF  number  of  wave  number  groups 

NALPHA  number  of  members  of  ALPHA  array 
NX1  number  of  profile  parameters  to  be  used 

in  zone  I 

NX2  number  of  profile  parameters  to  be  used  in 

zone  I  and  zone  II 
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5)  NAM4 

WNS  array  of  largest  wave  number  of  wave  number 

group 

6)  NAM5 

IF  index  array  of  wave  number  group 

7)  NAM  6 

EX1  array  of  profile  parameter  used  for  zone  I 

8)  NAM7 

EX2  array  of  profile  parameter  used  for  zone  I 

and  zone  II 

9)  NAM  8 

CTEMP  array  of  temperature  0°K  to  20000°K  in 
1000°K  increments 

10)  NAM9 

SIGMA  array  of  electrical  conductivity  data 
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B.4  FORTRAN  LISTING,  ARCRAD  III  AND  SUBROUTINES 


-j 


01/03/94 


CC»*0H  ALPHA!  A)  .HtQflEU29},BlirrTI64}«Cr2Tt20)  ,cos*»  t  e »  P 
1&3  MKO.f  *T<  >«CP<?U,  >i)75i£t  7kJ0)  aft(40i4Q.i}  »  S  l  A.  4  f>  »  , 

2TEHPR4C  2U  J  tWWB120»7fr>iRl«i7HjTt2J).rglu"i(?i>*CT6«PtZHtS10mT|?0> 
•P  ,  y  C  HN  t  « 1  ,t  tJD  .►‘+TlhT|S0i,S#C«*,21  J  i  *HI«:  .7  86. 4  J  yv*0<20  3  . 
•aiGtZO.ACPlMXCjO.ZAKQKlOJ  .O*ET{20,  1<J),  OEM  (T120.1CJ.  0x8512#.  10)  * 
■  ~*eUSU"(2D)  .SONEtUO  >  .SQ£h  tT(2g>  ,$CA«»(  20)  .  IF{lO)  .EX1I1CKEX2I19)  , 
*£<{101 .aN5llO}.gPAQI10) »Ci»iO|lfl1 
CliHMCV  AC i A  J  f  *  2  4  A A 
CCmmO*  /LIBEL/  NT  .Mol  2B  ) ,  TE"P(20> 

1  roppif ijHi ) 

2  FORMAT  IiJ|4A  I  _ _ _ 

1  rCRMiTHl) 

V  rQftr  a  T  C>XSMRAP.«*4HTEWP.7MSW50ENT,B)t5wS0AgS.  B X5 MSB*,  j  % 9  * 4 Mwc.CN ,  JO 
lKl44€«l5l*NP«TI*tn 

to  format  (1m  3U»ll*iT»lX>) 

14  r  ORM  a  T  J  1h  |i>ElK:}.2l’i**i.?I,A41PEU,i<2X>) 

13  rQRMATtiW03xi*«J£"P.l4x4MWvNCI  ■  11  *9*1*  XE*  OR) 

1*  rORnAt- iiHOlOxie^^JARfeNT  .  iP£14  7 ,  J.0X/HE  ABB  ■  1P€]4<71 

_ 17  rORMAT  IlMPlQltJHSIBSUM  a  !Pfl4.7> _ 

U  FORMAT  ClHgiOxi2«USWIAT)Ow  •  1PE1*,7) 

_ 13  FORMAT  t'  3,0^14,7^14,  7#ilA,7,El*,7j 

2V  FORMAT  1JN  4*eftAD«4xa(Eli,4'3xl> 

21  rORMAT  Uh09X44TkMPil4x?HCVl  _ 

22  FORMAT  JIM  2  C 1P|14 ,7,11)1 

21  FORMAT  11m  DX9h«N  1  fBt0,l23qPE11.4A2F)  | 

24  rORMAT  JlNOJ^AMTfeMp.lJxeMtTASuH) 

29  TOrmaT  1 1  h  JtiMt- ,  7,  1X71 

24  FORMAT  (1M  2MMNf  AMJF  6,0»9X)  I 

i ?  rOR^RT  i;h  BmCQR * ^ .A* B « t U 1 4 , z " H 

ll  FORmIT  |  liHOlDAPhSURAO  a  IREl*  7  .  lnx9NSC0HAu  a  1PE1*,7> 

_ NAMfc^  1ST  TA,MA>RA>TNALL>RNi  TCL/N4M2/Ty,0ALRMA|F.  IP«,  jIGTA, 

THZ/NAMl/N  I  r  ,N ALPHA ,NX1  ,Nx2/NaM4  /«  N  5/ N  AM*/ ]  f/  N  A M4  t  \  %  1?  N  A  *7  /  f  it 
2/Nam&/CTEnR/Nam*/MGma/Nam1Q/aI).A1*a2 
|f  ( 5ETF 1 T (  2  ) )  70*  724.7# 

70  DO  71  WJHajiM 

71  RE  AO ( 9 , 1?  )  )M  CNXXl  |F(NX>,U  ,F[NxX,2)  ,T  |NXX,3|  ,F  1*11,4) 

_  *£  40  l9.lLNt,43! _ 

71  READ  ti.2l  |F  NT (LI .Lai .12) 

SH*  SFaMD«xFI,T> _ 

SfcaO  (9«naMi> 

PEaO  I5,»AM2I _ 

me aq  iiiKTiiii 


M  4  - 


RE A 0  19,****) 
READ  I9.NAH7J 
PEAO  I3.MAMB) 

pCaO  ISAaMj 

RE  AO  I9.NAM10I 

pI I TE  U.fciun 

FRITE  li.NAMji 

pine  <#,m**J) 

JrflTEtiiHAPA ) 

-»nnF  T#.(*i«3) 

pR|T6  I».«ahA) 


2 


17 

* 

H 

4 


08/03/6# 

NaI*<  -  fcF  N  sOlRCE  STaTEPe^t  -  lf*|S)  - 

51i7f  e~~ iiT«*AM/  j 

NUf’F  ( A  , NAMfl ) 

pSltn*.*A«fll 

MBUE  U  .  N  A  M  t  f  I  > 

l  1 2 1  •  !  )  3  •  1 

A3a*0*M"lO,  *1.2*100  • 

PO  1  MU 

N.»  73 _ _  _ 

N*800*O 

I f  [RPfcAO  CU,BOOR6C,NMBC0)>75.T4#/5 

7 a  c*ll  Error  (j,i*hpqf  rm)li  PEaO  «ec» 

79  B  a  3 1  ^ a 0 

If  SRRbAOB  (|#.UN'N-Si«im.74i7« 

/»  nm8IN“U _ 

If  lRREAOai  1#.#KC*HPB  110)77. 7«*  7  7 
7  f  iHRf|]NaC 

|T  CR4tAUB(l4,REC*PMBiNH  tl,  74,76 
70  Ciy,.  CLOSE  U* ) 

C‘LL  RcaV  u*7 

)P$Ta4 _ 

TE"P|i)a271. 

RiQtl>*»A.R^ 

T  T • 1 90 C  , 

_ -D0-ifl-Wta.-iJ2i _ 

TT,TT*lyoc , 

TfKP£N)aTTa300,  _ _ 

Bfl~ R AO  t  N 1  **T* *1 1  *  C  {  T  T*Ta  )  /  £TMALL-r*  )  )*i  RR*1  .  >1 
RL*RaO)2I/HA 
TL*?900, 

CCN?Ta*»»S  l  <?.«^KJ_4)/l>t.^TWALL>  )»KRR.Rl>/IPm*WLM 
1-  UTL/tU-T-ALLM*iDL^twN.RL))»ALOBCTL/TRALt|-l  .)) 

DC  79  1*1,79 

- s  WiT." 1 7 -xa N?T*AkCt  |PST  >  - - 

79  |PSTa[R|1*» 

SCalp^*D*Lfm*.i, 7433j926-2 
NiNGLEaNALPMA*! 

ANOLEaO, 

sl*<f<Jil5Vx _ , _ — _ _ _ _ - 

00  ifQ  aA  ■  i  ,  *1 A  hSLfc 
AMlLE*AN«Lt*ROALRM 
$  1  N A  C  NAlafttm  AiaOLtl 
1»U  OStNA  ;  Nl-1  MS  f  MA  (  NA  J-SI  RA  NA-J  ) 

ALPHA  E  lPV“ *0A|>H/2, 

COSNAfl)  .  CO^lALRWAU)) 

DO  200  NA  a  2  j NaLPmA 

alpha  IP* 1>*lPH*X»*»1)»R0alPw 
sou  CO 5 N *T)» *>«COS[ALP^A)hA>I 

[ZO^EaD _ _  .  _ _ _  „ 

lNa[|2 

IM«IA«1 

- nITTTTi -  - - - 

b:»M2  _ 

px-JiXl 

JIO  ?05  { X  a  1 ,  N X 
ED  I)J)a*XlUX) 

210  DO  71*  |Jltli*K _ _  _ _ 


41 
4? 
41 
4  * 


50 

51 

91 

•  0 

44 

47 

4« 


00 


44 

102 

no 
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B.4  FORTRAN  LISTING,  ARCRAD  III  AND  SUBROUTINES  (Cont'd) 


H4  |  N 


fcru  S0u»cfc  Statist  -  ITNCS* 


0VO3'*4 


09/53/4* 

*»I1  •  tf*  90J«C§  ST4TE*t«HT  -  irstSl  - 


00 


T*4iT4/luOQ 


C  i  riM^  ITM.UZ.  I  I«ltTCIiL*E>f  ut.n*» 

||T?.2*NT 

NTjf.s  I-VK 

*T(,fcT.l 

■T?1«1T?-UC _ 

DQ  24 1  "«  | H , N  r 

r|  (  TfcHP(S)**MO, 4  9fl«€  M2  **-2931993 

«*i 

PO  320  NM.Nf2 
DO  320  HJ.l.NT? 
fC  3?0  ha.I.sa^hA 
TT^T  JIlfc.ftJiiAitO 

DO  UP  «■["!. <41 
H( *s-l 
JO  >J0 

im  r*(no  5  st  .LOS**  <  N*  > 


_a5_li0  _J8  _ _ 

350  4(N.NJ.NA|B.4p»T*(*H<T,a2.*(«AOIN.I|*a2-4*0tNl*<2j!*a.3 

_ &C_  400  *l«4Ti,N!22 

■  S»nT2-«i*i 

*U**W 

qC  400  «*.l  ,N*Lf"* 

0C  U1  Hj.lh.HS 

9]M.*|Ht-*2M«<LMSJ)'<2-H*0[*.S>*M< 
i  I N  t  •  0  j  N  .: 

_  B<P|,*fJ,NAMA|*T*9^T 

3?*  COM  IMJfe 

_  ir  ifjsi-srii  sao , *oo ,  *oc 

)i»  DD  J90  Hj.S5t.NT 

pt*Ta4  t«l*aZ*(*A(H*lJ)*a*-»40<«t)*l*2l 
fr  teiNTI  403. 3*9. 395 
_JB»  IJ^faS. Nl.*,4 
395  »<  N.Hj,  *A  i  ■*  t*T -D  1ST 

_  . 

U< N,kj2,Ni(a4|HT-«(NT 
393  qost |sut 
403  C  OS  I ISul 

00  410  S*|N,*T 
0C  *10  I  ■  ( •  «#H 

_ aULMUWiiJBli _ _ 

kUftCOD 

[ f  !RWt*U(19.BC04tC .SN|LDj |  4j3.4il.4L9 
*11  C*tL  EmkuI  (J.lBwfcOf  *N(LF  PS;  *  0  *f  L  ( 

413  00  420  I  ■  1 1  N|J 
H*9  jMII 

( f  {R»t*UB<14.RfC.N4BlN) ) 4 20. 4 1*,420 
A|4  t: 4 l L  E*4d»  n.i*siBT~SHTtE  flt*0  *»Pl) 


<20  C*LL  CL 0 4 E  r 1 4 1 

1 1  i  |IOHk.n  *l3.«l2. 

412  jAm*  _ 

Ofi  TO  429 

413  |M4a(l2 _ 


133 


50  1*4 


/I  l*Jf 


if  204 


22* 

2?« 


234 

231 

2i« 


42O0  449  MISN.ST 

474  9S«JN*0_^ _ _ _ _ 

iP'fflRfclCItlft.AICiNUf |«|)  431.*3?.431 

_ 432  C*LL  i»«04  IJjl9Mt0f  “M^f  4ia0  *ECJ 

431  if  U9HTf«»CH  >-«fcC<H  (-1  .!-*»  433.433.43U 
4HLJE*U  clo*e  f  16  ? 

00  TO  «29 

433  M»«  J SaO _ 

If  (fl*fc4D6U4.fi|C  *  N-9  ISM  4  34.432.  434 

434  JPIT a  IFOI 

DO  4  Jf  1 4  j  .79 

I&OMJUn-rftftCUKtll _ 

4  3  T  lfST.IF9T.4 
<**  C4LL  lili 

00  449  Rain. ST 


|0  441  •tfftN.Rf 
_ lSrl 

444  »aC2T|Sl*aN(  HI 

— _ liLXy^ZL  .431.4491449 

449  (W«(H«1 


44f"f?  Ui#-1>  4J0, 4449.430 

4<*?  I *c 2 TJJ  >±h*lL11 _ 

<1.2 

_«_!«  *n  _ 


490  VaC2TlN|*tt**||0-ll 

_ 11*1" _ 

491  »2M«* 

_ X4iT2*iI _ 

V9a  tf 4*Vf 

.  Bi’l.-  lU+WJJL" 


1*4  7^14.  .  J 
00  440  1*11.74 
VC2TI N|*4N( | | 

_ _ \f  t  VI.  1  4903.4903.492 

<909  02. Va ^ 

M.V|»U 


21XL-33J33333-V/9,. 


lllVOVl 

*N«t , -.193M97- | V*3)*f. 31933333*0/9. 
tU.OL-O" 

SL*  OH 

eO  T0^ifr" 

492  D<UH.P. 

.Had 


*Hafl . 

|0  4?4 
4  *  *  4  S  •  1  , 


1 JK»1.4 


• *2/90. • V* /J 94 O..04/2 72140. I 


_ 494  DSu"aOSU".t (El<M  -*«■¥( )/Aaa*4|*M ( Ak.V»3. ( *AH«V*9  .  (a*Na*«4.  , 

ONaOSUHa 

_ _ EU*lk-fl* _ 

OLa^S 

499  MOflM,  |*1 |a|TA«Tfc"FftW(*l 
|T4lU*l"HafeTA5g«<*t*El4 
\f  (0-19,1  445.441,491 

440  cOHT|Nu« 

441  COhTIHuQ 


290 

293 

299 

242 


27* 


>19 
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B.4  FORTRAN  LISTING,  ARCRAD  III  AND  SUBROUTINES  (Cont'd) 


o 

vO 


*1  m  -  trh 


-5Qu*CE  ITiTiUtm 


im$i 


0«/fl37<>6 


•^l*  •  Er«  90u»rE  ir*TttttNT  -  ir^tS) 


DO  <15  l*liNu 

_  DC  463  WlM.Nf _ 

46J  t*0RH<N,  1(4*. 

DO  41)  *•1,111 _ 

_ 14? 

oa  <*?  non. *t 

_ B  P__J  - 1  t*H _ 

OtMITIft,  1R)*0 

_ _ C"l*  IWl-C  • 

4*7  ONfT'N, 

IQ  <4*  L"ftiiM£ _ 

COffiOt I*  < • 0 . 

<M  a?iDtlffi»ti. _ _  _ 

DO  9355  N»( N # 1 T 

_ _ _ 

N5*llT2*N"l 

_ 

00  9359  N**I  •NH.PN* 

_ iL-lRlXLthLiMll _ LLLl?i&j31> 

»IQ  *r.i(’2-K» 

_ ll?  DO  9399  Il'liNl' _ _ _ 

RH|aif I ( H-l l«i 
.-  BNf  *  tf 1 1* ) 

|Q  9315 

_ iam*wBi  ln,  n«^AD(Ni  _ 

E«t^»»3N*>r  ,NM*»ri5pTlu  I  { 

lr  LN/-l!_m^l4iJli _ _ 

Mi  (KoniniNf  i  (*bKaR 

917  C*LL  LOQ*  U«flfli£"tC€»«| 

f  ?***5lO.-tN 

_ CTBMUMD-CM _ x _ * _ _ _ _ 

9179  Q  E*f  |  T  (  w «  |l)lGEM|1(N<  )  *  I  *  tf|*  H  »Q  J  1  *t M  Mi  J  •  Th  A  NS  i 

_ UL_l|il^LL^iex520*»?5 _ _ 

111  TfUNS»l*AN51 

_ _80  ?a  U?  _ _ _ _ 

*20  TKQfl.D* 

_ BO  93*  NJal.m _ 

«a«W.NJ 

m*PQ'(NQ, 

ET»C0«»  TKfim*SHC« 

C*U  tOQH  <  TKOW,  TK,CTM 

- call  LOOK  Tll.dMH.ctM 

IHiWSjTf-m _ _ _ _ 

“ CTH6S6CT*-CfcP 

92*  Ol|5  I  NO,  |H)»34lll(NQ.|«l-y|M«DilNMNH*(7W*N*t-TR4Nll 
ill*  TiAHSiifJuNI 

9  39  aliOUHHQMQl  |  ■]  •NMe»«0S  I  H  A  I  NM  On  1*5 

cClU&t  1*1  .CeRiCI  |*l<»kf3*«*D*  INaThO-CTHOS 


vtiBR'*e* inS-i, ji**> 

_ I  Tf 51 «g 

fc*f Cp* CnFTKS*!,  Hi  1»‘>R0P  ( n S -  1  , 1  ) 

call  lco*  «iina,i-M,c£»o 
thassi  •! , 

at  Ml  M**-’  ,  t*>  «QfeM]  T  |NS*1 .  |R>.m9«4)*QSImA(NA|«T*an$1 
T  N  (JR  ■  C  . 

_ I JiElH _ 

■  Ni- 1 

>35  |r(ft|*  *»,**) j  7  »)} 

*(3/  irtNF-MTi  >44.^44  5>? 

)«4  KI,S12-N> 

GO  |C  90 

_ nest-i  _ 

Nr»*T2-*> 

M*  T4  0N>*^DN<|N(N.Nr,H4>-ll(P<|]Nr<NA))*r’RaM|Ki,I4 
ETH0N"l<Cs-F4pn 
C»LL  L»3C"  <  < M (3 R *  |W,LTK| 

- —CALi.  LQQ*  tU_KSj[t>t*  C^wi 

T»AM»  IM-FM 
C***N5  ■  CT*?-CtH 

5<8  0*8SIKIr  i(liut4MM  lKf"MB^R*CSlVAvHA'*[TSAMS(“T*4MSI 
TUlNSl tTWANS 
JP  4  mr .  |  m  n  ij2.6gu,»E> 

>>>i  |  a r » 

»f  *sr»  i 

|f  C  |TE3T-j|  762 

400  QfiAQI  |K  »»3*AIH  I  R  1*4  *8*4(1  S  J  Mi  {  si  |  «TN  A  H$ 

JTf  6CQ1  CCNaOI  |M  laLQNArJ|  J  «)  •  *08R»DS  1  Ml  LA  {"CrPANS 

DC  6005  N ■ I L  «  4  T 
scnm<m)»d. 

$CF  m  j  t ( 6  3  *0  . 

®0  09  SC«ff5tNMD, 

QC  60 1  N»lN*niT 
00  60 1  1**2, Njr 

ONfl IN. J«(>OcM||(N,  H»l-JApS<N,  ffi I 
SC  ME  ftN|  •^asttlAM  *  QnG  T  t  N-  |R> 

SC|n|Tcn|,SgtH(f  (>*(  .QtMJ  r<h,  J»| 
tai  saA?&{6i  •'JOAhiqm  ♦  j»gs{ h,  ju , 

SQBl  [} 

sca»«o*a. 

00  6Q1>  1  l • 2 , n I r 
SOaAO-SawiU.uoiCt i R I 
SCisiT^TCg**  «LGMA0i£R) 
qQ  60J  m»{6,«dT 
00  603 

If|B|G<M|-fc<jRMil«N(  j))  602,  ROJ*8^) 

BlgtHiatnaRMi j\, j I 
mvmO<NI  i«m  |  i 
<0J  CONMNUfc 

I  8,11 
M^ITf  >415» 

QO  634  M>(N,NT 

ta<  ir8|T£l  4  jgi  rtnP<M  ,»VNQ|M{  ,fiju<h> 

•RJTF  |*.?4) 


3*5 

3*7 


OQl> 


_ pa  6001  N»N T I .NT21 

liS«N  <  2-N*  1 

_ DO  60*1  **•!,  — L»*i 

10  6QV1  lN«*.h  I r 

_ *Nt.ir(i*.ii»i _ 

*WF.|f<I*l 

_ |C  »>fl  f4*Mt«NNf _ 


Of/03/66 


436 


496 

451 


9?7 

92* 

931 

>39 
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B.4  FORTRAN  LISTING,  ARCRAD  III  AND  SUBROUTINES  (Cord'd) 


>1*  |N 


ef*  BOuftcC  lT*T|NSN?  *  tfNltl 


0*/DV*» 


"M1*  *  feF*  «Oo*CS  IT ATfH^NT  -  JfN(|| 


G#/DJ?*+ 


<3 

O 


00  *3*9  N«j*i.NT 

•  0*3  H»|ft  ll>?)l  TSHP  <N>  it  T>SUW)N} 
if  <i“-i7  hj,&oi.*iq 
4U*  H|TE  i  •  *  1  I 
*«|TE  (•*?!■ 

oc  »o*  n • i . i j a i 

OV<N]'(iaNkTlNn/<l*lU<*'l«*2-4lO[N«l|«*?l"J»1419*2*) 

_ *0*  WITI  <♦.?/>  <§Na<N»  tPV<  N) 

«#> NO l»7?-l 
DU«N0ir:i 
t (N0|V>«(« 
f  J  |G-*(  *»0|  VI*$|BT* 


J»l 

_ dc  ro«  b»iia.«<T _  _ 

TTTT^TTnUm 

f  Ll*»  SiOWAlMQ  JV) 

VlNC«(*4U(N)^«i0(N*l)|/D| V 

N | 

00  «••  l*2.NJlv 

MiM**  |Ng  _ 

Hp-MaiT  CL  Mil'll. 

1*4  J«J*l 


30  IQ  All 

7iT-tCNT.|iLflBCHlJ  l-*L3|tCTtP*»(  J-|  »  I)  /  (  *  LUO  S  C  TE  M  <  J  |  I  -  A  t0tt  <  C  TENI>  { J*  J  i 

in _ 

rrrTNrnTrrn^rTXi.Tj^iTrwiTJi/it*«*<  j-it  >)»^CNT**LObi5tow»i  j-ii  i 
n  i««r 


- BtrrBrrar 

41*  FS  IQ"!  t  I  HSU**  ( 

Hi  CTStlvue 


_ L^*Ai _ _ _ _  _  _ 

Do  •  »*  1 1  ■  1 «  nn 

4«>  SLN*svM»«„*rsiON*< r« i  i|»2,*rii0HM^-|ui, 

llOlMl«N$«v|NCM»^|fiM*(tl*SljM.4.afSl4iMA(ND|V-ll«FllaMi(NOlV(|/J< 

E1g,,*G  . 

hoS>*j. 
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B.4  FORTRAN  LISTING,  ARCRAD  III  AND  SUBROUTINES  (Cont'd) 
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B.4  FORTRAN  LISTING,  ARCRAD  III  AND  SUBROUTINES  (Cont'd) 


-  fcFN  5QUHCC  JTATfNfWT  •  IMlOl 


01/03/** 


SUHOVTIH  tOOii 

tnm*  «kM0f»i7icfl«fci sinforn  (Mi  ih  } itiMium  * - 

?TEM*iiii}t  *«»i  ?»i  rit*>mi<Ti?ii  »ff  i«i**  i  ?i  r*CfiH#ili>*lilfiitiii> 

**ICi79T7l*wRiT>&»7Vi  i«»  .o«Tiioao» *olnjMio*ioi  .silsili.ilti 
^ii5u^i^rUNHj/j>,^MTiao]ao*onai.  iFiHKtii  mi  ft*iU-U4— 

ro«o* 

I*  <*-.U  ?5.ZQ.> 

t.  1UHi:HS^3.H 

;>  (»  i  *«•»»•<»■ _ _ _  _ 

20  1 1 N C « . C 2 


23  5 ■ IX-Kl  (NX)  I/AINC 

‘  «*r<  Hi  •U..»Z*r{NKla.fi4|-].1>nNH'31«9«<S-t.>*lfi-.»»«Ft"l**4) 

*  ■  *  — - -  L  ■ 

i*  fuuf*\*^*i*fn*f*  _ _ _ _ _ _ _ _ 

00  TO  2* 

JLillil?!*! _  _ _ 

HER** 

*lNC«-DI _ _ _ _  _ 

_ _ SO  IP  » _  _ 

•0  A  1  *#C  •  .1 

_ *H»x/*i»c+i*i _  _ _ _ 

$0  TO  Z» 

»0  ilNTfl, _ _ 

Ai*f/il»c*  J4. 

i«i¥-»xt  iahm/ainc  __  _____ 

r  **i  ifTTt**,  i  i  •  cs^Tf  <**,  ?  >  a*  <  j*i .  i« i *- .1 i » 

■  i  h;m;5f - ~ 

33  r t. i  ao* ii«x*o«y>»zi  *f  * 

- ^  “  ~ J  - 

r t.4 j«»  i _ _ _ _ 

3t  flETUS* 

?o  r«*j, 

OETiM*  ’ 


12 


AEDC-TR-66-258 


APPENDIX  C 
ANNULUS  BOUNDARIES 
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C.l  WEIGHTED  AVERAGE  TEMPERATURE 


In  order  to  avoid  interpolation  of  absorption  coefficients  for  tem¬ 
perature,  the  boundaries  of  each  annulus  were  so  calculated  that  the 

o 

weighted  average  annulus  temperature  was  an  exact  multiple  of  1000  K. 

The  FINITE  ANNULUS  computer  program  to  accomplish  this  was  an  im¬ 
proved  version  of  the  ANALYTIC  TEMPERATURE  PROFILE  program 
described  in  Appendix  B  of  Ref.  Ma.  1*  FINITE  ANNULUS  was  developed 
in  its  present  form  on  the  small  computer  but  eventually  incorporated  into 
the  main  program,  Even  though  the  arc  as  a  whole  is  not  optically  thin,  the 
assumption  that  radiation  from  an  individual  annulus  varies  as  the  product 
of  annulus  volume  and  is  reasonable.  Thus  the  weighted  average  tem¬ 
perature  of  annulus  i  can  be  defined  as 


T.  = 

i 


(C.l) 


y 

The  exponent  on  T  was  determined  by  comparing  curves  of  T  for  values  of 
y  from  4  to  10  with  curves  of  pressure -dependent  total  radiance  which  had 
been  cross -plotted  from  density-temperature  results  of  Nardone,  et.  al. 

(Na.  1).  No  single  exponent  matches  all  the  curves,  of  course,  but  the  best 
fit  is  in  the  range  from  7  to  9  and  numerical  calculations  showed  that  results 
are  insensitive  to  choice  of  y,  Table  C.l,  The  exponent  was  therefore 
arbitrarily  fixed  at  8P  The  calculation  is  started  at  the  outer  boundary  (i  =  1) 
where  temperature  and  radius  are  presumed  known,  and  proceeds  inward  to 
the  centerline. 


TABLE  C.l 

COMPARISON  OF  CALCULATED  ANNULUS  RADII  FOR  SEVERAL 
VALUES  OF  TEMPERATURE  EXPONENT  (x  -  4,  rG=  1  cm,  Tc  -  14,000°K) 


T(°K) 


z 

5000 

7000 

9000 

11000 

13000 

14000 

7 

.9293 

.  8847 

.  8200 

.  7451 

.  6305 

.3805 

8 

.9311 

.8847 

.8212 

.  7454 

.6319 

.  3822 

9 

.9325 

.  8850 

.  8222 

.  7458 

.  6330 

.  3842 

%  Variation 

.3% 

.05% 

.3% 

.  1% 

.4% 

1% 
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Since  the  desired  values  of  T.  are  predetermined,  the  unknown  is  the 
inner  radius  of  the  annulus,  .  Annulus  radius  appears  both  explicitly 

and  as  a  limit  of  integration,  and  must  therefore  be  found  by  iteration  as 
described  in  Ref.  Ma.  1.  The  procedure  involves  selection  of  two  trial  values, 
interpolation  (or  extrapolation)  to  a  third  trial  value  and,  if  necessary, 
repetition  of  the  interpolation  step  using  results  of  the  latest  two  trials  until 
matches  its  preselected  value  within  _+  2°K. 


In  FINITE  ANNULUS  an  exception  was  made  for  the  core  and  inner¬ 
most  annulus.  Equation  C(  1  was  modified  to 


(TN)y*rN 


)  + 


(C.2) 


where  subscript  N  refers  to  the  core.  For  this  special  case  the  core  radius, 
only  appears  explicitly  and  can  be  solved  for  without  iteration. 

A  minimum  annulus  thickness  was  also  established*  If  a  calculated 
annulus  thickness  was  less  than  5%  or  ra,  that  average  temperature  was 
automatically  omitted  and  the  annulus  volume  included  in  the  next  higher  tem¬ 
perature  step,  A  listing  of  the  FINITE  ANNULUS  program  is  included  in 
Appendix  B  as  subroutine  FINANN  of  ARCRAD  III. 

C.2  OUTERMOST  ANNULUS 

o 

Below  3000  K  radiant  emission  from  the  arc  is  totally  negligible. 
However,  absorption  of  vacuum  ultra  violet  radiation  is  still  very  important. 
This  absorption  is  proportional  to  particle  density,  so  it  is  convenient  to 
compute  an  average  particle  density  for  the  outside  annulus  in  Zone  II,  where 
the  temperature  ranges  from  2500  K  at  the  inner  boundary  to  the  assumed 
wall  temperature  of  500  K, 

By  as  sumption,  the  temperature  profile  in  this  region  is  linear. 


Subscripts  L  and  W  refer  to  values  at  the  inner  and  outer  radii,  respectively. 
The  average  particle  density  is  given  by 


N  = 


v 


Lo 

V 


I  (p/po)  dV 
■'V 


(C.  4) 
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+  19  3 

where  Lo  is  the  Loschmidt  Number  (2.  69  x  10  particles /cm  at  STP), 
Assuming  that  air  at  these  temperatures  is  a  perfect  gas,  and  making  use 
of  the  cylindrical  geometry 


N  = 


(C.  5) 


With  equation  C.  3  this  integrates  to 


(C.6) 


Absorption  coefficients  are  then  directly  calculated  from  values  tabulated  for 
density  corresponding  to  STP 


k 


(k) 


STP 


X. 


N 

Lo 


(C.7) 


where  X.  is  the  mol  fraction  of  species  i  (0_,  N  or  NO). 

1.  O  Cj 
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APPENDIX  D 

TOTAL  RADIATION  CALORIMETER  MODIFICATION  TO 
RADIANT  INTERCHANGE  EQUATIONS 

The  sensing  element  of  the  calorimeter  designed  for  measurements 
of  total  radiation  was  ring  shaped  to  maintain  rotational  symmetry*  However, 
in  order  to  minimize  the  effect  of  thermal  conduction,  the  ring  was  recessed 
and  therefore  exposed  only  to  a  relatively  short  length  of  the  arc. 


Applying  the  radiant  interchange  integral,  Eq.  (2.  15),  to  this  situation, 
integration  in  the  4>  direction  is  unchanged  but  now  both  the  upper  and  lower 
limits  of  the  0  integral  are  functions  of  axial  position,  z,  and  an  axial  integra¬ 
tion  is  necessary.  Fig.  44. 


$ 

x  'n,m 


=  (Zirr  )(<JT  ) 


D(n)  -  D(n  +  An)  sin^  +  ~  sin ^ 


_2 

ir 


.  2 

sin 


©a© 


-(kR)t  / sin  © 
e-  +  e 


[(k  R )  +  <kR)  ] /sin© 


(D.  1) 


Here  h  is  the  half  width  of  the  calorimeter  ring  and  (z/h)  is  the  greatest 
relative  axial  displacement  at  which  any  part  of  the  calorimeter  sensing 
surface  is  visible.  Definition  of  half -width  and  choice  of  negative  (z/h)  were 
made  for  convenience  in  working  with  the  angle  variables.  The  former  changes 
the  symmetry  factor  on  the  integral  from  4  in  Eq.  (2.  15)  to  2  in  Eq.  (D.  1), 

The  dependence  of  limits  a  and  $  on  (z/h)  is  illustrated  in  Fig.  44* 

Two  cases  are  distinguished,  one  where  the  shield  overlaps  the  sensing  ring 
and  one  where  it  does  not.  In  both  cases  the  axial  integration  is  divided  into 
3  zones  and,  for  each  zone,  the  angle  limits  are  given  in  Table  D*  1  in  terms 
of  their  tangents. 

As  before,  the  integral  can  be  written  as  the  sum  of  two  integrals 
identical  in  form 


f 

c 


fO  *0  (z/h) 

max 


*/s“9Sin29ded<z/h> 


(D.  2) 
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The  required  double  numerical  integration  was  carried  out  on  a  digital 
computer.  For  the  geometries  of  interest,  the  value  of  the  integral  turns 
out  to  be  nearly  independent  of  the  location  of  the  radiating  surface  relative 
to  constrictor  wall  (dimension  b.  Fig.  44).  Therefore  the  ratio  of  the 
calorimeter  integral  to  the  complete  integral,  Eq.  (2.  16),  is  plotted  in 
Fig.  45  as  a  function  of  absorption  length  only.  A  polynomial  was  fitted  to 
this  curve  so  that  radiation  reaching  the  calorimeter  could  be  computed  as 
an  extra  step  in  the  ARCRAD  III  program.  The  ratio  of  the  summation  over 
all  annuli  and  all  frequencies  of  radiation  reaching  the  calorimeter  to  the 
equivalent  summation  of  radiation  leaving  the  arc  is  the  term  C  in  Eq.  (4.  1). 
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TABLE  D-l 


Note:  -z  is  a  positive  number  since  integration  is  from  -(z/h)  to  zero. 

rri3i  x 


a.  Geometry 


AIOIA875 


b.  Typical  Absorption  Path 
Fig*  1  Absorption  by  an  Infinitely  Long  Half-Cylinder 
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-  EMISSION,  INCLUDING  SELF-ABSORPTION 

- TRANSMISSION  -  ABSORPTION 


a.  Side  View,  Section 


b.  End  View 


Fig*  2  Typical  Emission  and  Absorption  Paths  to  be  Accounted  for  by  Radiant  Interchange  Model 
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ig.  3 


Examples  of  Subscripting  Nc 
Boundaries,  Annulus  Proper 
Sector  Radii  (N-  3,  m  =  3) 


k/c  (GM/CM. SEC.) 
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Fig.  6  Electrical  Conductivity  at  P  =  100  Atmospheres 


WAVENUMBER  (cm-1) 


AEDC-TR-66-258 


DATA  SOURCE 

I  I  NARDONE  ET  AL. 
r - SHERMAN  AND  KULANDER 

^ ^  PLUS  FREE-FREE  tv, 

I //  A  SCHULTZ  ET  AL.  KBT 


.2<V<  15  INCLUDES  MORE 
THAN  99%  OF  BLACK 
BODY  EMISSION 


Fig.  7  Range  and  Sources  (or  Absorption  Coefficients 
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ABSORPTION  COEFFICIENT  <CM~> ) 


0-10,000‘K 
x-(50OO 
■  -2O0OO 


b.  59,000-200/000  cm"1 


Fig.  8  Continuum  Absorption  Coefficient  for  Air  ot  100  Atmospheres  ond  Temperatures  of  lO/GOO^K,  15,OOOnK  and  20,000°K 
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CUMULATIVE  SPECTRAL  RADIANCE  (  WATTS  /  CM  3  -  STER  I 
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Fig.  11  Solf-Absorption  in  a  Uniform  Temperature  Arc,  ra  =  .1  cm  T  =  17,000° K 
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Fig.  12  Experimental  Setup.  The  High  Pressure  Air  Arc  is  an  the  Left.  On  the 

Right  is  the  Optical  System  Leading  to  the  Spectrograph,  a  Small  Part  of 
which  con  be  seen,  Lower  Right  of  Center 


Fig.  13  Optical  System  with  Rotating  Mirrar  (Rectangle  in  Center  af  Large 
Plexiglas  Disc)  ond  Associated  Electronics 
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Fig.  14  Schematic  Representation  of  Experimental  Apporatus,  Showing  Arc  Chamber,  Spectrograph,  Drum  Camera, 
Associated  Optical  Equipment,  Power  Supply  and  Synchronization  Equipment 
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Fig.  15  Design  Curve  for  Copper  Constrictor  Showing  Time  to  Surface  Melting  Temperature  as  a  Function 
of  Wall  Heat  Flux 
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Fig.  16  Exploded  Views  of  a  Copper  Cascade  Constrictor  and  an  Ablation  Type  Plexiglas  Constrictor 
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METHOD  OF  ASSEMBLY 


ELECTRODE 
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Fig.  17  Schematic  View  of  Copper  Constrictor  and  Components 
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Copper  Disc  Before  Etching 


Copper  Disc  After  Etching 


Fig.  18 


Copper  Disc  of  a  Cascade  Constrictor  Before  and  After  Etching.  Observe  thot 
the  Loose  Copper  Particles  and  Machine  Marks  are  Eliminated  by  Etching. 
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Copper  Disc  After  a  Symmetric  Operation 


Copper  Disc  After  an  Unsymmetric  Operation 


Fig*  19  Copper  Discs  of  a  Cascade  Type  Constrictor  After  Symmetric  and 
Unsymmetric  Heating  by  Arc 
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Fig.  20  High  Pressure  Arc  Chamber  and  Plexiglas  Holder  with  Copper  Constrictor 
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Fig.  21  Schematic  View  of  High  Pressure  Arc  Chamber  Assembly 
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Fig.  22  Detoils  of  High  Pressure  Apparatus: 

a)  Side  View  Showing  Current  Path  Through  Arc  and  Return  Through 
High  Pressure  Vessel;  b)  Top  View  Showing  Opticol  Path  ond 
Apertures  Added  to  Reduce  Stray  Light 
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1  WEAKLY  DAMPED  LC  CIRCUIT 

2  RC  CIRCUIT 

3  STRONGLY  DAMPED  LC  CIRCUIT 


Fig*  23  The  Power  Pulse  and  Its  Components 
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Fig-  25  Arrangement  of  Apertures  in  the  Spectrograph 
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Fig.  26  Intensity  of  the  Carbon  Standard 
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Fig.  27  Arc  Cage  with  Calorimeter  Installed  and  Calorimeter,  Partially  Disassembled 
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©  COPPER  RING  SENSING  ELEMENT 
©  ALUMINA  CERAMIC 
©  PLEXIGLAS  SPACER  DISC 
©  COPPER  DISC 
©  THERMOCOUPLE 


NOTE:  INSIDE  SURFACE  OF  SENSING 
ELEMENT  THREADED  ON 
MODEL  3  ONLY 


Fig.  28  Section  Drawing  of  Total  Radiation  Colorimeter 
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Fig.  29  Colorimeter  Response  to  Constant  Heat  Flux  of  5  KW/cm2  for 
0-71  ms,  Followed  by  Adiabatic  Isolation 


105 


INTENSITY  (ERG/CM -A-STER) 


TEMPERATURE  (°K) 


Fig*  30  Intensity  of  Continuous  Radiotion  of  High  Pressure  Air  Plasma  in  the  Temperature  Range  12*2O,OO0°K* 
Pressure  100  Atm;  Wave  Length  8330A;  Assumed  Thickness  of  Radiating  Layer  1/10  mm 
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Fig.  31  Typical  Instabilities  of  a  High  Pressure  Arc  Column.  Constrictor 

Diameter:  5  mm;  Pressure:  31  atm;  Gas:  Nitrogen;  Current:  700  Amp 
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Fig.  33  Discrete  and  “Smear”  Pictures  of  a  400  Amp,  100  Atm,  Air  Arc  Showing  Fluctuations  of  Arc  in  Constrictor 
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Fig.  34  Drum  Camero  Film  of  o  100  Atm  Air  Arc  of  a  Current  of  115  Amp  Showing  Formation  of  an  Arc  Cor# 
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Core  is  Farmed 


Fig.  35  Magnified  Sections  from  Drum  Camera  Film  Fig.  34.  Characteristic 
Phases  af  Arc  Development  are  Shown  in  Detail 
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Expansion  of  Arc  Column 


Diffuse  Arc  Column 


Care  is  Formed 


Fig.  36  Discrete  Pictures,  Taken  Simultaneous  with  ''Smear*'  Picture  of  the  115 
Amp  100  Atm  Air  Arc,  Fig.  34.  Kink  Type  Instabilities  Grow  Stranger 
With  Time 
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Fig.  37 


High  Speed  Pictures  of  an  Ablation  Type  Arc,  Showing  Its  Remarkable  Stability*  Constrictor;  Plexiglas; 
Inner  Diameter:  3  mm;  Current:  400  Amp;  Pressure  150  Atm 
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CURVES:  H  Q  =  H (9500  °K) 


Fig.  38  Air  Are  Centerline  Temperoture  Inferred  from  Electrical  Conductivity 
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TEMPERATURE  (°K) 
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Fig.  39  Measured  Temperature  Profile:  Current;  115  amp.  Pressure,  103 
Atmospheres;  Constrictor  Diameter  2.5  mm. 
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8330  CONTINIUM 


WAVE  LENGTH  A  7770  8450 

INFRA  RED  SPECTRUM  OF  HIGH 
PRESSURE  AIR  ARC 


SMEAR  PICTURE  OF  ARC  COLUMN 
DURING  SPECTROGRAPHIC  EXPOSURE 


Fig.  40  Spectrum  and  Drum  Picture  of  115  Amp  Arc.  Temperature  Profile  was 
Measured  Over  Width  Marked  A  in  Smear  Picture 
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Fig.  41  Curve  Fitted  to  Radial  Temperature  Distribution  of  a  High  pressure 
Air  Arc,  Fig.  39 
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CUMULATIVE  FRACTION  OF  TOTAL  RADIATION  LEAVING  ARC  COLUMN 
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Fig.  42  Cumulative  Net  Radiation  Based  on  Measured  Temperature  Profile, 
p  =  100  atm,  I  =  100  amp,  TCL  =  17,700°K 
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Fig.  43  Colorimeter  Traces;  P  *  100  Atm,  1  *  410  Amperes 
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o.  Shield  Does  Not  Overlap  Sensing  Element  (S  >  0) 


ZONE  I 


b.  Shield  Overlaps  Sensing  Element  (S  <  0) 
Fig.  44  Calorimeter  View  of  Arc  Column 
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Fi^.  45  Ratio  of  Calorimeter  Integral  to  Complete  Transmission  Integral 
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